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Preface 


Thermodynamics deals with energy conversions and predicts the directions of 
these energy transformations. Thermodynamics is needed to understand many 
processes on Earth, be they physical, chemical, or biological. Thermodynamics 
is critical to study the atmosphere (lapse rate, foehn winds, circulation), hydro- 
sphere (latent and sensible heat, pressure dependence of freezing/boiling points), 
geosphere (geothermal gradients, mineral stability), and the biosphere (redox 
zonation, evolution of biogeochemical cycles). 

Thermodynamics is based on two major laws. The first law concerns the conserva- 
tion of energy and the second law states that energy is increasingly being dispersed. 
These two principal laws of thermodynamics are general and simple, both in words 
and equations, but their applications may be challenging. Citing the Feynman lectures 
on Physics (Feynman et al., 1963): “Thermodynamics is a rather difficult and complex 
subject when we come to apply it, and it is not appropriate for us to go very far into 
the applications in this course. ...The subject of thermodynamics is complicated 
because there are so many different ways of describing the same thing” It is also 
challenging to teach to undergraduate students, because some might better grasp it 
with a rigorous approach involving calculus and derivation of each equation, while 
others would understand the topic better using everyday examples and limiting the 
use of mathematics. 

This introduction to thermodynamics and equilibria is part of the first-year chem- 
istry class for Earth (System) Science students at Utrecht University. It aims to prepare 
students for more advanced second-year classes in physical chemistry, mineralogy, 
and petrology and to provide the basic concepts of relevance for atmospheric, marine, 
climate, and environmental sciences. Mathematics is limited to the inevitable and 
notations are kept as simple as possible to maintain accessibility of the text. This 
introduction is limited to ideal, equilibrium systems, neglects the temperature and 
pressure dependence of enthalpy, entropy, and heat capacities, and does not explic- 
itly mention the reversible or irreversible nature, or the path, of system changes. 
These and other simplifications at the expense of scientific rigor and comprehensive- 
ness are hopefully compensated by students’ improved understanding how to use 
thermodynamics for problems in Earth System Science. 


vi Preface 


This concise, accessible introduction was written during a one-month workation 
in the Ardennes (May 2023) with the aim to have a revised first-year course ready 
for the 2024 spring semester. The author acknowledges Jeannot Trampert and in 
particular Lubos Polerecky for critical feedback on the first draft, teaching assistants 
Chris Bil and Thomas Sanders for checking the second draft, and Ton Markus for 
upgrading the quality of the figures. 
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Part I 
Thermodynamics 


Chapter 1 A) 
Introduction ERE 


Abstract This chapter introduces the macroscopic approach to matter and defines 
thermodynamic systems and their exchange with the surroundings. The state vari- 
ables temperature and pressure are presented and linked via the ideal gas law (an 
equation of state for gas). 


Keywords Partial pressure - Temperature - Ideal gas law - Thermodynamic system 


Thermodynamics is a scientific discipline firmly rooted in physics, chemistry and 
engineering with implications for the universe and thus biology, earth sciences 
and astronomy. It describes the behavior of matter on a macroscale, sets the rules 
for conversions of energy from one form into another, and governs the direction 
of processes. Engineers and physicists applying thermodynamics focus on energy 
conversions and conservation to develop efficient engines, refrigerators, or heat 
pumps. Chemists use thermodynamics to study chemical reactions and phase equi- 
libria. Thermodynamics is key to understand many processes on Earth including 
adiabatic compression of air (Fohn winds), the lapse rate due to expansion of air, 
geothermal gradients in the inner Earth, the effect of pressure on freezing and boiling 
points, the strength of atmospheric circulation, the stability of minerals as a function 
of pressure, temperature, and composition, redox zonation in marine sediments and 
the evolution of biogeochemical cycles. 

Thermodynamics is a science about the macroscopic world, i.e., the human scale 
and larger. Most of the thermodynamic concepts have been developed, and can be 
understood, without knowledge of the molecular structure. However, sometimes it is 
easier understood if we use our knowledge on the properties of individual molecules. 

Thermodynamics is based upon a few statements, i.e., thermodynamic laws, that 
are based on years of observations and multiple clever experiments connected through 
theory. Most chemical thermodynamics textbooks present four laws, while engi- 
neering treatments are often limited to the first and second laws. The first law is the 
conservation of energy: energy can be converted from one form into another, but 
overall, no energy is lost or gained. The second law is concerned with spontaneous 
change and states that energy is, overall and over time, increasingly being dispersed. 
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4 1 Introduction 


Alternatively formulated, the second law reads that spontaneous processes are those 
which increase the entropy of the universe. The term entropy will be introduced and 
explained later. 


1.1 Macroscopic Approach, Yet Everything is Made 
of Atoms 


Chemists deal with atoms, ions, and molecules and aim to explain the properties 
of matter and predict the synthesis of compounds using atom-level knowledge. 
Although thermodynamics is a core subject in academic chemistry curricula, it deals 
with macroscopic systems, i.e., large amounts of matter rather than a few molecules. 
Matter exists in three phases: solid, liquid or gas that differ in the number of parti- 
cles (atoms, molecules, ions) per volume, hence the distance and interactions among 
particles, and their spatial orientations and distribution (Fig. 1.1). In a gas molecules 
are far apart and move randomly through a volume that is primarily empty space. 
A gas is thus homogenous, fluid, and compressible. Moreover, for an ideal gas the 
identity and interaction of the particles can be ignored at the macroscopic level. In 
a liquid the particles are close to each other but are free to move relative to each 
other, thus a liquid is fluid but less compressible. In a solid constituent, atoms, 
ions, and molecules are generally close to each other and in fixed positions. More- 
over, their identity governs the physical and chemical properties of the solid. Pure 
substances, matter with a homogenous and definite chemical composition, may exist 
as a solid, liquid or gas depending on the external conditions (temperature, pres- 
sure). For instance, at the Earth surface, H2O, water, may exist as solid (ice), liquid 
(water) or gas (water vapor). A single phase may contain multiple substances, e.g., 
air consists of nitrogen, oxygen, carbon dioxide, etc. 


States of matter 
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Solid Liquid Gas 


Fig. 1.1 The most common states of matter: solid, liquid and gas 
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-= ae Ln 
1 bar 10° Pa 
1 mm Hg 133.32 Pa 
1 torricille or torr 133.32 Pa 
1 atm 101.325 Pa 


Thermodynamic theory is general and applies to solid, liquid, and gas phases. 
However, most thermodynamics has been developed based on the behavior of gas. 
Consider a container filled with a gas. Gas exerts a pressure on the container walls, 
and this can be understood in terms of randomly moving particles that collide with the 
wall. The result of these collisions is a force perpendicular to the wall. Gas pressure 
(P) is thus defined as: 


pe (1.1) 
~ A’ . 


where F = force (N) and A = area (m°). Pressure acts equally in all directions in a 
gas and is thus a scalar quantity. Pressure has the unit Nm~? or Pa (the SI unit), but 
multiple other units are commonly used (Table 1.1). 

Because of the large amount of space between the molecules in a gas (Fig. 1.1), 
and thus limited interactions, gases with different compositions will mix well. If they 
do not react, the total pressure of the gas mixture is the sum of the partial pressure 
of each gas contributing to the mixture. The partial pressure of gas A (Pa) is the 
pressure it would have if it were alone. Dalton’s law states that the total pressure 
(Piot) for a mixture of gases A, B, C, is given by: 


Prot = Pa + Pg + Pc (1.2) 


For instance, for air (78% N2, 21% O2) with a total pressure of 1 atm, the partial 
pressure of nitrogen and oxygen gases would be 0.78 and 0.21 atm, respectively. 
The particles moving freely in a gas have kinetic energy. Temperature is a measure 
of how much kinetic energy each particle has on average. The higher the tempera- 
ture, the more energy a system has, all other factors being the same. Temperature 
is not a form of energy, but a measurable parameter to compare amounts of energy 
of different systems. Most solids, liquids, and gases expand roughly linearly with 
increasing temperature because temperature governs the kinetic energy of particles. 
The expansion of liquids such as mercury and alcohol are commonly used to quan- 
tify temperature using the freezing (0) and boiling (100) points of water for calibra- 
tion. This centigrade scale has been superseded by the closely related Celsius scale 
anchored to the triple point of water (0.01 °C, see Sect. 4.3) at which ice, water and 
gaseous water are at equilibrium. The Celsius scale is widely used but the scale is 
arbitrary: it is based on water and its macroscopic phases. The thermodynamic, or 
absolute, temperature scale is independent of the substance used in the thermometer 
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and goes down to the minimum possible temperature of —273.15 °C. The absolute 
scale, or Kelvin scale, is expressed in K (without a degree sign). Temperature in 
degrees Celsius and Kelvins are related: 


K = °C 4597345 (1.3) 


This Kelvin scale should be used in thermodynamics. 

Having defined pressure (P) and temperature (T). We can define two reference 
states: standard temperature and pressure (STP) refers to P = 1 bar and T = 273.15 
K = 0.0 °C, while standard ambient temperature and pressure (SATP) refers to T = 
298.15 K = 25 °C and 1 bar for P (1 bar = 10° Pa). 


1.2 Ideal Gas Law 


Experimental studies on the physical properties of gases have resulted in a few 
empirical laws that were eventually combined into what we call today the ideal gas 
law. Boyle (1627—1691) studied the relationship between the pressure and volume 
of a fixed amount of gas at constant temperature and observed that the product of 
pressure and volume is constant. Boyle’s law states that the volume (V) of a given 
amount (e.g., mass) of a gas is inversely proportional to pressure (P) when the 
temperature is constant: 


V x —, (1.4a) 


where « represents the proportionality symbol, or alternatively: 
P . V = constant. (1.4b) 
Subsequent work by Charles (1746—1823) showed that at constant pressure (P) 
the volume (V) of a given amount of a gas is directly proportional to its absolute 
temperature (T). Charles’s law states: 


VaT. (1.5a) 


Or alternatively, 
= = constant. (1.5b) 


Next, Avogrado (1776-1856) related volume (V) and the amount (n in moles), 
but at fixed temperature and pressure. Avogrado’s law states: 
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Van, (1.6a) 


Or alternatively, 


V 
— = constant. (1.6b) 
n 


These three gas laws can be combined because volume (V) appears in all three: 


pa (1.7) 
YS is 
P 


The proportionality symbol (x) can be turned into an equality through the 
introduction of a proportionality constant (R): 


yap (1.8a) 
= R. —, 8a 
P 
This equation is usually written as: 
PV =nRT, (1.8b) 


which is the well-known ideal gas law. R is the universal gas law constant, which 
has the value of 8.314 J mol“! K~! in SI units, when n is expressed in moles. Note 
that thermodynamic data are often provided in various units (e.g., pressure in bar, 
atm or tor rather than Pa; volume in liters rather than m?); the units and values of R 
then need modification to maintain consistency of units. For instance, the value of R 
is 0.0820568 when expressed in L atm mol~! K7!. 


1.3 System, Surrounding and Equations of State 


Systems are a central concept in thermodynamics. The object of interest is defined as 
the system while everything else, the rest of the universe, is defined as surroundings. 
A system can be separated from the surroundings via imaginary boundaries or phys- 
ically real boundaries such as the wall of a container. The system can be a chemical 
reaction taking place in a solution (e.g., calcite dissolution), it can be your cup of hot 
coffee cooling while you are enthusiastically reading this text, or the entire System 
Earth receiving short-wave radiation from the sun and emitting long-wave radiation 
to the surrounding cosmos. 

A system interacts with its environment via mass transfer, or energy transfer via 
heat or work exchange (Fig. 1.2). An isolated system does not exchange heat, work, 
or matter with the surroundings. A closed system does not exchange matter but can 
exchange energy with the environment. An open system exchanges matter with the 
surroundings and may exchange energy (as heat or work) as well. Changes to the 
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system are called isothermal if the temperature is kept constant, isobaric if pressure 
is kept constant, isochoric or isovolumetric if the volume is constant, and adiabatic 
if no heat is exchanged, i.e., if it occurs in thermal isolation (Fig. 1.3). 

The state of a macroscopic system is described using observable quantities, i.e. 
state variables, such as pressure, temperature, volume, and number of moles. Equa- 
tions of state relate the various state variables of a system. The ideal gas law (Eqs. 1.8a 
and 1.8b) is an equation of state for an ideal gas. For one mol of a gas (n = 1), the 
system is fully determined with two out of the three variables (V, T or P). For instance, 
for a standard temperature and pressure (STP, 0 °C and 1 bar) the gas would then 


Open Closed Isolated 


wA A 
> = 


Fig. 1.2 Open systems exchange matter and energy; Closed systems exchange energy but no matter; 
Isolated system can exchange energy nor matter 
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Pressure (atm) 


T= 300 


isochoric 
1 1 

1 1,25 1,5 1,75 2 2,25 2,5 
Volume (L) 


Fig. 1.3 Isothermal, adiabatic, isochoric and isobaric changes to a thermodynamic system (e.g., a 
gas) 
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have a volume of 22.4 L. For SATP conditions (25 °C and 1 bar) the molar volume 
(Vm) of an ideal gas would be about 24.5 L mol7!. 

Thermodynamic variables are referred to as intensive if they are independent of 
the amount, while extensive variables are additive, i.e., they depend on the sample 
size. Of the four variables in the gas law (Eqs. 1.8a and 1.8b), P and T are independent 
of the amount of gas and are referred to as intensive variables, while V and n are 
proportional to the amount of gas and are extensive variables. In thermodynamics, 
changes in extensive quantities are associated with changes in the respective specific 
intensive quantities and their product has the dimension of energy: e.g., PV in the 
ideal gas law. The ratio of two extensive variables can be an intensive variable, e.g., 
density (mass divided by volume). 

The energy of a system is related to all the other measurables of the system 
via equations of state. Thermodynamics literally means ‘heat movement’ because 
it describes how the energy of a system relates to measurable variables. If the state 
of a system shows no tendency to change, i.e., the system is at equilibrium, the 
state functions have values which are independent of the history of the system. 
Consequently, the changes in a function of state do not depend on the route by which 
one goes from one state to another, i.e., they are path independent. Equilibrium 
thermodynamics, the topic of this course, focuses on differences in system states and 
we use the symbol A to indicate changes in state variables and functions. 


Box 1 Math intermezzo: Partial derivatives and state functions 
Equations of state usually involve multiple variables. The total differential of 
a function F (x,y,z) is defined as: 


OF OF OF 
GUE = || =] a@bese || = ]}) @hy se (| SS dz (1.9) 
ax yz dy X,Z dz x,y 


where the derivative of F is taken with respect to one variable at a time while the 
aF 


others are kept constant. The first term (+) ae is the derivative of F with respect 
tox only andis the partial derivative. These partial derivatives sometimes reveal 
relationships between state variables and are often used to formally define basic 
properties such as heat capacity, compressibility, etc. 

As an example, suppose we aim to quantify the pressure dependence on 
temperature of an ideal gas, assuming that the volume and number of molecules 


remain constant. The relevant partial derivative is 


oP 
OT J vn 


We start with rewriting the ideal gas law (Eqs. 1.8a and 1.8b) to isolate 
pressure on one side of the equation: 
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nRT 


Next, we take the derivative of both sides with respect to temperature T, 
while considering n and V constant, and obtain 


oP ð (nRT nR ð nR 
E = Se pnn, (1.10) 
ƏT Jynn TNV V oT V 


In this way, we have analytically derived how pressure varies with temper- 
ature and have revealed the phenomenological law of Guy-Lussac stating that, 
for a given mass (n) and volume (V), the pressure is proportional to absolute 
temperature, i.e., AP/AT = constant. 
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Chapter 2 A) 
The First Law: Work, Heat creck ie 
and Thermochemistry 


Abstract This chapter presents the first law of thermodynamics in terms of heat 
exchange and work. Internal energy, heat capacity and enthalpy are defined and 
enthalpy changes during reactions and phase changes are discussed. Latent and 
sensible heat are introduced, and adiabatic processes are related to the atmospheric 
lapse rate, potential temperature, and the geothermal gradient. 


Keywords Zeroth Law - First Law - Heat - Work - Enthalpy - Heat capacity - 
Latent heat + Sensible heat - Adiabatic processes + Phase changes 


Before elaborating on the first law and its implications, it is necessary to introduce 
the zeroth law of thermodynamics and to further define our systems. Consider two 
systems, A and B, that are closed, i.e., no matter can exchange, but with different 
temperatures. If these are brought in contact, heat will flow spontaneously from the 
hotter to the colder system till temperatures become equal, i.e., thermal equilibrium is 
reached. Energy transfer from one system to another due to temperature differences 
is called heat. If a third closed system, C, is in thermal equilibrium with system A, 
then Tc =T; and system C must be in thermal equilibrium with system B also. This 
illustrates the zeroth law of thermodynamics stating that 


if two systems are in thermal equilibrium with each other and a third system is in thermal 
equilibrium with one of them, then it is in thermal equilibrium with the other also. 


Systems exchange matter and/or energy with their surroundings. Thermody- 
namics has adopted a system-centric view, i.e., matter added to, heat delivered to, 
and work done on the system have a positive sign. Conversely, work done by the 
system or heat flowing out of the system to the surroundings have a negative sign 
(Fig. 2.1). 
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Fig. 2.1 A closed system ° 

interacts with its S u rro u n d | ngs 
environment via work (w) 
and heat (q) exchange. Heat 
added to and work done on 
the system are positive by 


+w -w 
m A a 
System 
Surroundings 
2.1 Work 


Energy is defined as the capacity to do work. Work is being done when motion occurs 
upon the action of a force. Work, w, is quantified as the product of force, F, and the 
distance, d, over which it acts: 


w=F.d. (2.1) 


Work has the unit of Joule (1 J = 1 Nm), like other forms of energy. In thermody- 
namics a common form of work is expansion/compression work because of volume 
changes of the system. Consider a cylinder with a piston (Fig. 2.2). Any process that 
increases the gas volume within the cylinder under the piston will push the piston 
upwards and in this way the system performs work on the surroundings, i.e., it loses 
some energy in the form of work. The amount of work (w) done by the system for a 
change in volume (AV) is 


w = — Pex AV, (2.2) 


where Pex is the external pressure of the surroundings, and the minus sign is needed 
to satisfy the convention that work done by the system should be negative (Fig. 2.1). 
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1.4 atm 


a 
> “a 
h h AV 


| Viinal aa 
A 


Fig. 2.2 Expansion work. Left: A piston exerts 2.7 atm on the gas within the cylinder and the gas 
has volume (Vinjriai). Middle/right: Following lowering of the external pressure to 1.4 atm, the 
volume of the gas increases (AV = Vfinal — Vinitial = h(eight) x A(rea)) because of work done 
by the system w = —P.,, AV 


2.2 Heat 


Heat, also called thermal energy, is denoted by the letter q and has the unit of joule. 
Traditionally, the calorie (cal) unit has been used for heat. One calorie is the energy 
needed to heat 1 mL of water from 15 to 16 °C, which is equal to 4.184 J. The amount 
of heat needed to change the temperature of a system depends on the temperature 
change, AT, and the heat capacity (C): 


q=C-T (2.3) 


The heat capacity C is an extensive property that depends on the amount of 
material in the system. The molar heat capacity, Cm, is the amount of heat needed 
to raise the temperature of 1 mol by 1 K and has the unit J mol`! K~!, while the 
specific heat capacity, C;, is the amount of heat for 1 K temperature increase per 
gram of material with the unit J g7! K~!. Using the molar or specific heat capacities 
Equation (2.3) would then respectively become: 


q =nC,, AT, (2.4) 
and 

q =mC,AT (2.5) 
where m is mass (g) and n is the number of moles. As an example, consider the 


energy required to heat 18 gram (1 mole) of water with a specific heat capacity (Cs 
of 4.18 J g7! K7!) from 25 to 50 °C. Using Eq. (2.5), we then arrive at 1882 J. 
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2.3 The First Law 


The first law of thermodynamics can be formulated in multiple ways but in fact 
states: energy can neither be created nor destroyed, i.e., energy is conserved. The 
total energy of a system is defined as the internal energy for which we use the symbol 
U. In all phases (gas, liquid, solid) atoms and molecules are in motion. Internal 
energy comprises all types of energy: e.g., kinetic energy from molecular rotation 
and bond vibration, and potential energy from intermolecular attraction and chemical 
bonds. Because we cannot measure all types of energy in a system, the absolute total 
internal energy cannot be known. However, differences in internal energy (AU) can 
be measured and calculated. 

An isolated system does not exchange matter or energy with the surroundings; 
hence the total internal energy of the system does not change. This leads to the 
statement of the first law of thermodynamics: 

For an isolated system, the total energy of the system remains constant. 


Mathematically, this means: AU = 0 (2.6) 


This statement of the first law has limited utility because it is for an isolated 
system. An alternative, more useful way to express the first law (for a closed system) 
is the change in internal energy is the sum of heat (q) added and work (w) done on 
the system, mathematically: 


AU=q+w (2.7) 


This is the formulation of the first law used in this course. The internal energy U 
is a state function, i.e., it is path independent, while the heat g and work w are not 
state functions and depend on the path. In more advanced treatments, this equation 
is usually presented in calculus form as dU = ôq + dw, where d represents an 
infinitesimally small change and ô is used to articulate that infinitesimally small 
changes in q and w depend on the path. 

Another useful way of phrasing the first law is to consider expansion work 
(w=-PAYV): 


AU =q — PAV (2.8) 


If we add heat to a system with a constant volume, e.g., a container, there is no 
P-V work, AV=0, then 


AU =q+0=yqy, (2.9) 
where qy is the heat transfer at constant volume. Accordingly, the change in internal 


energy for a system at a constant volume is equal to the heat exchange that can 
be measured. This is very useful for engineers and others working with volume 
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constrained systems (containers, reactors). However, for earth scientists studying 
processes at a certain pressure, or chemists studying reactions in the laboratory at 
atmospheric pressure, it is more useful to quantify the change in heat at a constant 
pressure, i.e., the enthalpy change, AH: 


AH = qp = AU + PAV, (2.10) 


where qp is the heat transfer at constant pressure. This equation directly follows 
from Eq. (2.8) and leads to the formal definition of enthalpy, H, as 


H=U+PV (2.11) 


where the symbol = is used for a formal definition. 

Although the change in internal energy (AU) is the more fundamental quantity 
(used by physicists), the change in enthalpy (AH) is usually easier to measure and 
thus more commonly used. At constant P, the difference between AU and AH is 
due to the PAV term (Eq. 2.10). Volume changes involving liquid or solid phases 
are small, hence the PAV can be ignored and AH AU. However, volume changes 
can be large if gases are involved and then enthalpy changes are a little higher than 
internal energy changes (2.48 kJ per mol at 25 °C). Enthalpy changes follow the 
sign convention introduced earlier (Fig. 2.1): processes releasing heat are exothermic 
and have a negative AH because the system loses energy, heat uptake by a system 
corresponds to positive AH values and are called endothermic. 

The heat capacity (C) was defined as the amount of heat required to increase the 
temperature of a substance by 1K. Similar to heat (q), the heat capacity of a gas also 
depends on conditions, i.e., whether the process occurs under constant volume or 
constant pressure conditions. In case of constant volume, all the heat will increase 
the internal energy, and thus the temperature, because no expansion work is done. At 
a constant pressure, some of the heat goes into work and thus more energy is needed 
to cause the same temperature rise. The constant pressure heat capacity Cp is thus 
higher than the constant volume heat capacity Cy. Mathematically, Cy relates AT 
to the change in internal energy: 


AU = qy = CyAT (2.12) 

whereas Cp links AT to the enthalpy change 
AH = qp = CpAT. (2.13) 
Volume change with heating of liquids and solids is generally limited and their 
Cy and Cp are very close in value and any difference can usually be neglected. For 


an ideal gas, Mayer’s relation links Cy and Cp: 


Cp = Cy +nR (2.14) 
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Having discussed the internal energy changes of a gas at a constant temperature 
and constant pressure, we now consider an isothermal change. Isothermal conditions 
imply that AT is zero, hence AU = 0 (Eq. 2.12). The first law then reads 


AU =0=q+w 
and consequently, for isothermal processes 
q =—w (2.15) 


An isothermal compression (decrease in volume) causes a net outflow of heat, 
while isothermal expansion (increase in volume) must be compensated by heat uptake 
to maintain constant temperature. 


2.4 Thermochemistry: Enthalpy Changes in Chemical 
Reactions 


All chemical substances have enthalpy stored in the form of chemical bonds. When 
a chemical reaction occurs, the change in enthalpy, AH, is equal to the total enthalpy 
of products (final condition with new chemical bonds) minus the total enthalpy from 
the reactants (initial condition with pre-existing chemical bonds) 


A,H = A products — Hreactants, (2.16) 


where A,#H is the enthalpy change of a reaction, as indicated by the subscript 
r. However, absolute enthalpies of substance cannot be determined (like internal 
energy), only relative values, i.e., changes in enthalpy can be determined. In other 
words, we need a set of standards against which enthalpies of reaction can be 
measured. The standard enthalpy of formation of a certain compound, A H°, is 
defined as the enthalpy change for the formation of one mole of that compound from 
its elements in the most stable form under standard conditions (STAP, 25 °C and 1 
bar). The subscript ¢ is added to indicate that is the enthalpy of formation and the 
superscript ° that it relates to pure substances at SATP. For example, the gas N2 and 
the solid graphite are the most stable form of elements N and C at standard conditions 
and their enthalpies of formation (AH?) have been set to zero. Thermochemists have 
constructed consistent databases with standard enthalpies of formation that can be 
used to calculate the enthalpy of reaction at standard conditions and making use of 
reaction stoichiometry. 

Specifically, for the reaction of œ moles of substance A and f moles of substance 
B to form y moles of substance C and ô moles of substance D 


aA +BB —>yC+8D (2.17) 


the enthalpy of reaction at SATP (A, H°) is calculated as follows, 
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A-H? = [y x AH#(C) +8 x AH?(D)] — [a x AH?(A) + B x AH?(B)] 
(2.18) 


or more generalized: 
AH? = 5 vu; AH? (products) — 5 v; AH; (reactants), (2.19) 


where v; is the stoichiometric coefficient of substance i in the reaction. 

To illustrate this method to calculate A, H°, consider the oxidation of 1 mole of 
sucrose (C12H22011; s) to CO, (aq) and water (1) at 1 bar. The relevant reaction and 
enthalpy of formations are: 


Reaction | C12 M22011 (s) + 1202(g) => | 12CO2(aq)+ 11M,0(/) 
AHọ —2226.1 kJ mole! |0 kJ mole™! —412.9 kJ mole™! | —258.8 kJ mole7! 


and the enthalpy of reaction is calculated using 2.18: 


AH? = [12A H; (C O2(aq)) + 11A H? (O (D)] 
— [AH}(Ci2H2011(8)) + 12 x AH? (02(8))] 
= 12 x —412.9 + 11 x —285.8 — 2226.1 — 12 x 0 = —5873 kJ 


This enthalpy of reaction is also known as enthalpy of combustion or heat of 
combustion, i.e., the enthalpy of change when one mole of a compound reacts 
completely with excess oxygen gas. When using thermodynamic data, one should 
note that each compound and each phase of a certain compound have distinct internal 
energies and enthalpies and that the phase [solid (s), liquid (1), dissolved (aq) and 
gas (g)] should therefore explicitly be indicated to use the correct energy. This is 
because, as we will see below, phase transitions cause enthalpy changes. 

This is one way to obtain enthalpy of reactions, the other method is based on 
Hess’s law. 

Since enthalpy is a state function, only the initial and final conditions do matter, 
and the pathway (intermediate reactions) does not matter. Hess’s law states that the 
overall enthalpy change for a reaction is equal to the sum of enthalpy changes for 
the individual steps in the reaction. 

To illustrate Hess’s law, consider the Haber(-Bosch) process, the common way to 
produce ammonia (needed for fertilizers to produce food for the world population) 
in large quantities from nitrogen and hydrogen gas. It is a multiple-step process: 


Step 1 : 2H2(g) + N2(g) > N2H4(8) AH, =? 
Step 2 : No H4(g) + Ho(g) > 2NA3(g) AH = —187.6kJ 
Combined : 3H3(g) + No(g) > 2NA3(g) AAj42 = —92.2kJ 


(net reaction) 
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Hess’s law implies that AH, + AH; = AHAj+2; hence AH, = AHi — AA = 
+95.4kJ. The first step is thus endothermic consuming heat (the system/reaction 
receives heat), while the overall reaction is exothermic and produces heat. 

The heat of combustion or enthalpy of combustion is usually presented in kJ g~! 
(i.e. specific combustion enthalpy) rather than kJ mol~', the unit used for enthalpy 


of formation (anp). This specific combustion enthalpy is a measure of the energy 
density. 

While the heats of combustion of methane, ethane, propane, butane and pentane 
are —890, —1560, —2220, —2877 and —3509 kJ mol`! and differ substantially, 
their respective specific combustion enthalpies are more alike: —55.5, —51.9, —50.4, 
—49.5 and —48.6 kJ g~!, respectively. For comparison, typical specific combustion 
enthalpies in kJ g7! of kerosene (—46.2), diesel (—44.2), crude oil (—43) are also in 
this range, while those of anthracite coal (—32.5), lignite coal (—15), wood (—15 to 
—20) and peat (—15 to —20) are much lower. 


2.5 Enthalpy Changes During Phase Changes 


Matter exists in three phases (solid, liquid and gas) and transformation from one 
phase to another involves heat transfer, i.e., a change in enthalpy (Fig. 2.3). The tran- 
sition from solid to liquid, melting or fusion, requires addition of heat (endothermic 
process). The enthalpy of fusion AH fusion is defined as the amount of heat neces- 
sary to melt a substance at its melting point without changing its temperature at 
constant pressure. Standard enthalpies of fusion AH7,,,,,, are available at 1 bar and 
25 °C. Freezing or solidification, the opposite process of melting, is exothermic and 
the enthalpy change has the same magnitude as that of fusion but the opposite sign 
(negative for freezing, positive for fusion). 

Vaporization (or boiling) is the transformation of a liquid into gas or vapor and the 
enthalpy of vaporization A Hyqp is the energy required at constant pressure to vaporize 
one mole of pure liquid at its boiling temperature. Condensation, the opposite of 
vaporization, has the same enthalpy as vaporization but the opposite sign (negative 
for condensation, positive for vaporization). 

Sublimation is the transformation of a solid into a gas. The enthalpy of sublimation 
A HĦ;ubı is the amount of heat required to convert a solid into a gas without going 
through the liquid phase. An example is dry ice (solid carbon dioxide) which at 
atmospheric pressure goes directly into a gas without formation of liquid carbon 
dioxide at a temperature of —78.5 °C. Condensation or deposition, the reverse process 
has again a different sign for the enthalpy change than sublimation. At a constant 
temperature, enthalpies of sublimation can be directly calculated from the enthalpies 
of fusion and vaporization using Hess’s law: 


A Fu pi = AF fusion F A Hvap (2.20) 
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Fig. 2.3 Phase transformations involving solid, liquid and gas. Melting, sublimation and evapora- 
tion involve heat consumption, while freezing, deposition and condensation release heat 
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2.6 Turning Ice into Steam: Latent and Sensible Heat 


Water plays a major role in System Earth: it dominates the cryosphere (and obviously 
the hydrosphere) and impacts the biosphere, geosphere and atmosphere. We have 
covered enough material to calculate the temperature response during the heating of 
1 gram of ice with an initial temperature of —20 °C to a steam of 120 °C at earth 
surface pressure of 1 atm (Fig. 2.4). 

Upon addition of heat to ice of —20 °C, the ice will heat up and with a specific 
heat capacity (C,) of 2.09 J g7! K~!. For one gram of ice, 41.8 J (20 x 2.09) will 
be required to reach the melting point of ice. At the melting point, 0 °C, further 
addition of heat will melt the ice but the system remains at the same temperature, 
and 331 J will be needed to convert one gram of ice into liquid (enthalpy of fusion). 
Further addition of heat will warm up the water and 418 J are required to reach 
the boiling point at 100 °C (since water has a specific heat capacity of 4.18 J g7! 
K~!). At the boiling point, liquid water will gradually turn into steam (evaporate), 
and this requires 2260 J (enthalpy of vaporization), but temperature does not change. 
Finally, further addition of heat will increase the temperature of the steam. To reach 
120 °C, this will require an additional 36.8 J, because steam has a specific heat 
content of 1.84 J g7! K~!. The total energy required is about 3088 J, of which 84 
% is used for phase changes (10.7 % for melting and 73.2 % for boiling) and does 
not result in a temperature increase. The energy consumed or released during phase 
transitions at a constant temperature is also known as ‘hidden’ or latent heat. This 
contrast with sensible heat, that results in a temperature change of the system. Latent 
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Fig. 2.4 Latent and sensible heat during warming of ice to steam 


heat fluxes are an important component of Earth’s surface energy budget, e.g., the 
evaporation/transpiration at the Earth surface and subsequent condensation of water 
in the troposphere. In the Earth’s interior, latent heat is also released if the liquid 
outer core crystallizes at the inner core boundary. 


2.7 Adiabatic Changes, Lapse Rate, Geothermal Gradient 
and Potential Temperature 


Adiabatic processes are thermodynamic state changes without heat transfer (q = 0). 
Accordingly, the first law AU = q + w then reduces to 


AU = w=-—PAV (2.21) 


stating that the change in internal energy of a system equals work. Adiabatic processes 
are usually associated with a temperature change of the system because of changes 
in volume. Compression (due to increasing pressure) causes temperature increases 
and expansion results in cooling due to reducing pressure. If a system expands or 
compresses very fast, or heat exchange occurs very slowly, it will not have time to 
exchange heat with its surroundings. Describing the system volume change as an 
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adiabatic process is then an appropriate approximation. For instance, air transport 
upward in the atmosphere can be modelled as an adiabatic process because the heat 
transfer between the air parcel and the rest of the atmosphere is slow on the timescale 
of its upward motion. A rising magma plume or air escaping after opening a bottle 
of champagne can also be treated as an adiabatic process. The geothermal gradient 
can be explained by adiabatic compression (Box 2). 

For an adiabatic process, the change in internal energy can be determined from 
two expressions (Eqs. 2.21 and 2.12): 


AU = —PAV and AU = CyAT 


In which Cy is the heat capacity at constant volume. 
Equating these two expressions 


CyAT = -PAV (2.22) 


shows that compression (AV < 0) leads to a positive AT, hence an increase in 
temperature, while expansion (AV > 0) results in a negative AT and is thus 
accompanied by a decrease in temperature. 

Adiabatic expansion can quantitatively explain the atmospheric lapse rate, i.e., 
the well-known decrease in temperature while moving upward in the atmosphere. 
Mathematically, lapse rate (T) is the temperature gradient: 


=e (2.23) 


where z is height (in m). 
To link the lapse rate with the first law, we treated an air parcel as a closed system 
and differentiate the ideal gas law for one mol (PV = RT): 


PAV + VAP = RAT 


and then re-arrange to isolate AT 


_ PAV+VAP 
z R 


AT 


Then substitute it into Eq. (2.22): 


PAV + VAP 
¢ (eee 


: )+Pav=o (2.24a) 


using Mayer’s relation (Eq. 2.14) R = Cp — Cy, we obtain 
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(= + VAP 


PAV=0 2.24b 
C Gy )+ (2.24b) 


Which can be rewritten as 


VAP 
Ce ` 


Cy 


PAV =— 


(2.25) 


This is the thermodynamic equation for adiabatic processes. It is sometimes 
Cc 


presented in integrated form as P V© = constant. Substituting Eq. (2.25) into the 
first law for adiabatic processes (Eq. 2.22) using the specific heat capacity at constant 
volume, we obtain 


VAP 
mCy AT = Cp 
Cy 
which can be rewritten as 
pCpAT = AP (2.26) 


where p is 7,, the density of the gas. 
Combining this equation with the equation for hydrostatic equilibrium of the 
atmosphere, 


AP =—pgAz (2.27) 


with g as the standard gravitation acceleration at the earth surface (9.8 m s7?) and z 
as height (m) to eliminate the pressure term, we eventually arrive at the lapse rate: 


since the specific heat capacity Cp value for dry air is about 1 kJ kg~! K~!. This 
theoretical estimate for the dry air is fully consistent with the experience-based lapse 
rate of 1 °C per 100 m (Fig. 2.5). 

However, air contains moisture that condenses upon cooling (cloud formation) 
and then releases (latent) heat (A H). This heat is added to the internal energy of the 
air. In this case we must extend Eq. (2.22) with latent heat: 


mCyAT =—PAV+AH (2.28) 
And since AH is positive (added to air), the adiabatic cooling (AT) will be 


less, thus the environmental lapse rate will be smaller, consistent with observations 
(Fig. 2.5). Moist adiabatic lapse rates are typically 0.5 °C per 100 m. 
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Fig. 2.5 Dry adiabatic lapse rate applies till level of condensation, above which latent heat is 
released (Eq. 2.28) and the lapse rate becomes less 


Pressure change induced adiabatic processes increase temperatures upon 
compression and decrease temperatures following expansion. Atmospheric and 
ocean scientists therefore use potential temperature, which is defined as the temper- 
ature that a parcel of air/water would attain if adiabatically brought to a standard 
reference pressure, i.e., earth surface conditions. Meteorologists use a standard pres- 
sure of 1000 mbar (1000 hPa) for reference, while oceanographers use the ocean 
surface pressure (0 dbar). Changes in temperature solely caused by compression or 
expansion are generally not of interest when studying atmosphere or ocean dynamics 
or heat content. 

The potential temperature (0) of an ideal gas is given by 


g=7( 7 g 2.29 
se ats 


where P is pressure, Po is the reference pressure, T is the current temperature in 
K, R is the gas constant and Cp is the specific heat capacity at constant pressure. 
Meteorologists use a value 0.286 for the ratio a Potential temperature (0) in the 
ocean is always lower than the actual temperature by about 0.1 K for every km 
depth increase. This limited increase of in situ water temperature relative to potential 


temperature is due to the large heat capacity of water. 
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Box 2 Geothermal gradient of the inner Earth 


The geothermal gradient, i.e., the temperature increase with depth into the 
Earth, can also be explained by adiabatic compression. The temperatures of 
the Earth’s mantle and outer core are close to the adiabatic temperature curve. 
For adiabatic compression, a similar analysis as for the atmospheric lapse rate 
due to expansion can adopted, starting with a modified form of Eq. (2.26) 


where m is mass, Cp is the specific heat capacity of the minerals, and œp is the 

fractional increase in volume per degree increase in temperature. Specifically, 

the volume change in atmospheric adiabatic processes a in Eq. (2.25) relates 
Cy 


to the volume change T Væp in the inner Earth adiabat (Eq. 2.30). 
Re-arranging this equation to isolate the adiabatic change of temperature 
with pressure: 


AT = TVap E Tap (2.31) 
AP E mCp z pCp l 
where p is the density 7. 


Using the appropriate hydrostatic pressure Eq. (2.27), we obtain after re- 
arrangement 


AT e gTap 
M Te 


(2.32) 


Using representative values for physical properties in the Earth’s interior 


(Table 2.1), we can calculate the adiabatic temperature gradients: at = 
0.88 Kkm™! at the core-mantle boundary and at = 0.29Kkm™! at the 


inner-core boundary. 


Table 2.1 Physical parameters in the outer and inner core near to the core-mantle boundary 
(CMB) and inner-core boundary 9ICB) (Lowry, 2011) 


Property Unit Core-mantle boundary | Inner core boundary 
Gravity, g ms 10.7 4.4 

Density, p kg m~? 9900 12,980 

Cp IK 'kg' |815 728 

T K 3700 5000 

Volume expansion |1076 Jkg7! 18.0 oi 

coefficient, œ p 
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Chapter 3 A) 
Entropy and the Second Law cieca; 


Abstract This chapter presents the second and third laws of thermodynamics, both 
at the macroscopic and microscopic scale. The state variable entropy is introduced, 
and it changes during reactions are quantified. A spontaneous process is recognized 
to cause an increase in entropy. 


Keyword Second Law - Third Law - Entropy - Spontaneous process 


The first law has shown to provide much insight into how nature functions, but it 
does not explain why some processes occur spontaneously. A spontaneous process 
proceeds in one direction, and once initiated, will continue. For instance, if you drop 
a rock from some table height, it will fall spontaneously, but a rock on the ground 
does not jump back to table height spontaneously. If we have two flasks, one at 
vacuum and the other one at 1 bar, and connect these, the pressure will equalize 
to 0.5 bar in each of the flasks, but without outside intervention gases will never 
unmix so that one flask will attain vacuum and the other 1 bar again. Heat will 
always flow from a warm to cold reservoir and not the other way around. These 
spontaneous processes cause a decrease in useable energy: gravitational potential 
energy for falling rock, work for the expanding gas and thermal energy for the heat 
flow. A process such as the dissolution of sodium chloride in water is spontaneous, 
yet it is an endothermic process (AH? = +3.88 kJ mol7!), consuming enthalpy 
and cooling the water. Apparently, changes in energy constrained by the first law 
are insufficient to explain the spontaneous dissolution of kitchen salt while you are 
cooking. Another thermodynamic law (the second) and a property called entropy are 
required. Entropy is a thermodynamic state function denoted by the symbol S, which 
is explained below. 
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3.1 The Second Law 


There are multiple formulations of the second law that emphasize its various implica- 
tions. The second law not only makes a statement of the direction of change, but also 
puts constraints on the efficiency of energy conservations. The latter aspects of the 
second law are often illustrated with a Carnot cycle of a heat engine and are primarily 
the domain of physicists and engineers, while chemists give more attention to the 
direction of reactions. While work can be converted to heat with 100% efficiency, this 
is not the case for the conversion of heat to work. There is an inherent loss of useable 
heat when it is used to perform work. This observation underlies Clausius’ statement 
of the 2nd law that heat can never pass from a colder to a warmer body without some 
other change, connected therewith, occurring at the same time. Another formulation 
for the 2nd law by Kelvin states there exists no process in which the sole result is the 
complete conversion of heat absorbed by a system into work. Clausius’ statement 
articulates the direction of spontaneous processes, while Kelvin’s statement empha- 
sizes the degradation of energy when it is transferred as heat (not all heat is available 
for work). While the first law states that the sum of heat and work (=internal energy) 
is conserved, the second law adds that the interconvertibility of heat and work is 
asymmetrical. Another, useable formulation of the second law states that sponta- 
neous processes are those which increase the entropy of the universe. Accordingly, 
we must properly define the entropy change. Entropy will first be illustrated using 
statistical thermodynamics (and on the way we will introduce the third law of ther- 
modynamics) and then the macroscopic entropy definition will be presented and 
elaborated. 


3.2 Microscopic View of Entropy and the Third Law 


The microscopic definition is based on statistical mechanics. Each macroscopic ther- 
modynamic state has a specific number of microstates, W, associated with it. A 
microstate is a particular way in which the total energy of the system is distributed 
among the individual molecules of the system. Consider a gas distributed over two 
bulbs, one is empty, the other one filled with four molecules (Fig. 3.1). After opening 
the valves between these bulbs, the molecules are re-distributed and 16 possible 
molecular arrangements can be found (=2+; 2 bulbs and 4 four molecules), each 
configuration is equally likely. The probability of the system returning to its original 
state is 1/16 (=2~*). Now consider the case of one mole of gas in one of the bulbs and 
the other bulb at vacuum again. One mole of gas contains 6.022 x 107 molecules, 
i.e., Avogrado’s number (Na). It will be clear that the probability for one mole of 
gas containing +1074 molecules to return to its original state in one of the bulbs only 
is extremely low (2-10), Accordingly, the number of molecular arrangements, or 
microstates W, increases during spontaneous changes. In other words, the random- 
ness or disorder of the system increases in spontaneous processes. Boltzmann defined 
the entropy as follows: 
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Fig. 3.1 Possible distribution of four molecules before and after connection of two gas bulbs 


S = kgln W. (3.1) 


This equation links entropy (S) with the Boltzmann constant (kg = 1.381 x 107% 
J K`!) and the number of microstates (W). The Boltzmann constant is via Avogrado’s 
constant (Na) related to the macroscopic gas constant (R): 


R=kg- Na 


This statistical thermodynamic view of entropy as a measure of randomness or 
disorder is instructive to explain some of its properties (e.g., its dependence on phase 
and temperature). Gas has more molecular randomness than liquids (because of the 
larger distance between molecules) and its entropy is therefore higher than that of 
liquids. Liquids allow movement of individual molecules and have therefore more 
randomness than solids, in which molecules are in fixed positions. Consequently, the 
entropy will decline when liquids freeze to form a solid phase. The dissolution of 
a well-organized sodium chloride crystal in water leads to increase in randomness 
because of the disruption of the crystal, but the hydration of the chloride and sodium 
ions somewhat lowers the randomness; the net result is that randomness and thus 
entropy increases, driving the dissolution of kitchen salt in water. 

Lowering of temperature decreases the kinetic energy of molecules and atoms. 
When the absolute temperature of a system approaches 0 K, then the motion of 
individuals molecules approaches zero and there is only one microstate (W = 1). 
Accordingly, application of Boltzmann’s law (Eq. 3.1) then indicates that the entropy 
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should be zero. This constraint is known as the third law of thermodynamics stating 
that the entropy of a perfect crystal at zero kelvin is zero. 


3.3 Macroscopic View of Entropy 


The macroscopic definition is based on a detailed analysis of heat-work cycles and 
the efficiency of heat engines (Carnot cycle). Carnot has shown that cyclic processes 
involving isothermal and adiabatic changes of an ideal gas cause no change in 
(internal) energy because of the return to initial conditions. Consider a heat engine 
with two reservoirs at temperature Tp (the source) and T, (the sink), receiving heat 
qn and transferring heat g,, respectively (Fig. 3.2). Because T, is larger than T, we 
can extract work (w) from the heat engine. 

Carnot showed that the maximum efficiency (e) of such a heat engine with fixed 
reservoir temperatures is 


Tn — T, 
e= = = ieee (3.2) 
dh Th 


The efficiency of an engine can thus only be 100% if T, is absolute zero. The 
Carnot, i.e., ideal, efficiency of an automobile engine having a Ta of 1100 K and an 


Fig. 3.2 A hot reservoir at 
Th supplies 20 kJ heat qn to a Hot source 
heat engine performing 5 kJ 
work (w) and the remaining 
15 kJ is lost as heat (qc) to 
the cold reservoir at Te. The 
efficiency is the work done 
divided by the heat supplied 
(qn) 
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ambient temperature T, of 295 K is about 73%; however, the average automobile 
engine is typically much less efficient due to other causes. 

Conservation of energy implies that 


w+ qnr + qe = 0. (3.3) 


Isolating w from this equation and inserting the expression in (3.2) then yields 


=——, (3.4a) 
dh Th 
or its equivalent 
dh dc 
—+—=0. 3.4b 
T, + T. (3.4b) 


Which implies that the ratio 4 is a state function whose change is zero in the 
Carnot cycle. 
Following Clausius, the entropy change can then be defined as 


AS = T (3.5) 


Next, we consider the spontaneous transfer of heat q from a hot reservoir (Th) 
to a cold reservoir (T}). The entropy decrease for the hot reservoir (4) is smaller 


than the entropy increase of the cold reservoir ( $ , and thus the overall entropy 


change AS is positive, as required by the second law (Fig. 3.3). Clausius inequality 
statement implies that for any spontaneous process 


AS > 0. (3.6) 


Fig. 3.3 Heat (q) 
spontaneously flows from a 
hot reservoir with 
temperature Ty to a cold 
reservoir at temperature Te. 
The entropy change AS is 
positive because the entropy 
decrease for the hot reservoir 


k is smaller than the 


ao 
AS= T 


entropy increase of the cold 


reservoir ( E ) : AS= 
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3.4 Entropy Changes During Reactions 


Since entropy is a state function, the change in entropy (A^, S) is the difference 
between the entropy of products and reactants: 


A,S = Sproducts — reactants 


Following the same approach as for enthalpy changes, we can obtain the entropy 
change (A, S?) for the reaction of œ moles of substance A and 8 moles of substance 
B to form y moles of substance C and ô moles of substance D 


aA+BB>yC+6D 
from the standard molar entropies at STAP ( 5°) as follows 
A,S° = [y x S°(C) +8 x S°(D)] — [a x S°(A) + B x AS°(B)] (3.7) 
or more generalized: 


A,S° = >, vi AS’ (products) — > v; AS° (reactants) (3.8) 


The standard molar entropy, S°, of a compound is the molar entropy of a 
substance in its standard state at 25 °C and 1 bar pressure. These values are tabulated 
at STAP conditions and normally expressed in J K~! mol`! (rather than kJ mol~! as 
enthalpy). Moreover, these are absolute values because the third law sets S° at zero 
at zero K. 

Consider the following reaction (Table 3.1): 


2H2(g) + O2(g) > 2M0 (1) (3.9) 
The entropy change of the system (=reaction) is calculated, 


A,S° = [2 x (m0 0)] — [2 x S°Hb(g)) + 1 x AS°(O2(g))] 
= —326.68 J K`! 


Table 3.1 Thermodynamic 

data for hydrogen and oxygen H2 (8) O2(g) H20 0) 

gas and liquid water AH}, kJ mol! 0 0 —285.83 
5°, J K7! mol! 130.68 205.14 69.91 
AG*, kJ mol"! 0 0 —237.13 
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There is a loss of entropy because three moles of gas (disordered) are turned into 
two moles of liquid water (less molecules and less disordered). Taking at face value, 
this negative entropy change of the system (A Ssys) suggests that water should not 
be stable and that the reverse reaction, transforming liquid water into hydrogen and 
oxygen gas, should be spontaneous. This is incorrect. The issue is that the second law 
states that spontaneous reactions should increase the total entropy of the universe, 
not that of the system only. To resolve this, we calculate the entropy change of the 
surroundings (A Ssurr ) realizing that the enthalpy change of the surrounding A H? ,, 
should be the negative of that of the system (A, H? = AHȘ,): 


AH urr = ~AHSys (3.10) 
Using the above thermodynamic data and Eq. (2.18), 
A, H? = 2(—285.83) — (2x 0 + 1 x 0) = — 571.66kKJ 


For constant P, we can use Eq. (3.5) (the entropy definition): 


AH?,, AHR, 571.66 x 1000 
AS), = 4 = Sour = 8 “= 1918.3 J K7! 
np iy 298 


The total entropy change (A Sror) is thus. 


AStot = AS$ 


0 + ASS, = —326.68 + 1918.3 = +1591.6 J K`! 


In conclusion, for a spontaneous process 


ASTor = A Ssys F ASsurr > 0, (3.11) 


Consistent with the 2nd law statement at the end of Sect. 3.1: spontaneous 
processes are those which increase the entropy of the universe. The entropy change 
of the system may be negative, but the total entropy change should be positive for 
any spontaneous process. This leads to a few simple statements: 


ASTot > Ofor a spontaneous process (3.12a) 
ASTot < 0 for a non - spontaneous process (3.12b) 
ASTot = 0 for an equilibrium (3.12c) 


Finally, any process taking place in the universe is spontaneous and leads to an 
increase of Sros. Consequently, the A Sros > 0 criterion implies a unique direction of 
time and underlies the maximum entropy production concept in use in Earth System 
and Ecosystem sciences. 
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Chapter 4 ®) 
The Gibbs Free Energy get 


Abstract This chapter introduces the Gibbs free energy and its relevance for the 
direction of phase changes and reactions. The dependence of Gibbs free energy on 
temperature and pressure is presented and used to derive the (Clausius-)Clapeyron 
relation. One-component phase diagrams and Gibb’s phase rule are presented. 


Keywords Gibbs free energy - Clapeyron relation - Gibbs’ phase rule - Phase 
diagram 


The second law provides a robust framework for predicting whether processes will 
be spontaneous or not. However, the assessment is made at the level of the universe 
rather than the system alone. To resolve this, the Gibbs free energy, another state 
function, must be introduced. 

Starting with Eq. 3.11: 


AST ot = ASsys + ASsurr >0 


we can replace A Ssurr using the definition of entropy (Eq. 3.5): 


ASsurr = Ter = ae 


Noting that the latter equality only applies under constant pressure conditions 
(see Sect. 2.3). The enthalpy change of the surroundings is the reverse of that of the 
system (Eq. 3.10): AH sy, = —AHs,;. Hence, we can obtain an equation for the 
system only: 


AH sys 


ASsys = >0 


which after re-arrangement reads 


AH sys = TAS sys <0 
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Table 4.1 Temperature and spontaneity of processes 


AH AS -TAS AG High T Low T 

= + = E Spontaneous Spontaneous 

= = + + Nonspontaneous Spontaneous 
+ a E Spontaneous Nonspontaneous 
— + + Nonspontaneous Nonspontaneous 


and the subscripts can be dropped to eventually obtain an equation for spontaneous 
reaction based on systems properties only 


AH — TAS <0. (4.1) 

The Gibbs (free) energy (G) is formally defined as 
G=H-TS (4.2) 
and is a state function because the constituent terms (H, T, S) are state variables. The 
change in Gibbs free energy (AG) can thus be evaluated based on the difference 


between final and initial states. 
Accordingly at constant temperature and pressure, 


AG = AH—-TAS. (4.3) 


This useful equation provides the basis for predicting spontaneity, hence the direc- 
tion of processes, a criterion for equilibrium and quantifies the energy available for 
work or metabolism of organisms. Comparing Eq. 4.3 with Eq. 3.12, it is clear that 
at constant P and T, 


AG < 0 for a spontaneous process (4.4a) 
AG > 0 for a non-spontaneous process (4.4b) 
AG = 0 for equilibrium, no tendency of the system to change (4.4c) 


The Gibbs free energy function also provides first-order guidance for the temper- 
ature dependence of chemical and phase equilibria. The change in Gibbs free energy 
AG = AH — TAS contains two terms: a change in enthalpy or heat of reaction 
(AH) and the product of temperature and the change in entropy (T AS). The temper- 
ature dependence of the direction of phase transformations and chemical reactions 
relate to this second term. Processes can be divided into four categories based on the 
temperature response (Table 4.1). 
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The temperature at which processes switch from being spontaneous to becoming 
non-spontaneous can be obtained from AG = AH — TAS = 0 which after re- 
arrangement leads to 


_ AH 


=e (4.5) 


The Gibbs function can also be interpreted in the context of efficiency analogue 
to potential energy in mechanics and the ratio of work done to heat supplied in 
energy analysis. The Gibbs free energy gives the maximal energy available to perform 
work, while the enthalpy gives the total amount of energy released due to phase 
transformations or chemical reactions and the T AS represents the loss term or 


non-useable energy needed for expansion work. The efficiency is then defined as 
AG 
AH: 
Similar to the introduction of standard enthalpies of formation (see Sect. 4.2), 


we can also define standard Gibbs free energies of formation (AG*“,) which is the 
change in Gibbs free energy for the formation of one mol of that compound from 
its elements under standard conditions (STAP, 25 °C and 1 bar). These AG‘ values 
are tabulated, reported in kJ mole! and provide information on the stability of 
compounds. For instance, the AG‘, values of water in liquid and gas forms are 
— 237.1 and — 228.6 kJ mol"! respectively, indicating that liquid water is the stable 
phase at 25 °C and 1 bar pressure. 

These AG‘, values can also be used to calculate the Gibbs free energy of a reaction 
(A,G?°) at 298 K and 1 bar: 


A,G? = > vu; AG* (products) — 5 v; AG‘ (reactants), (4.6) 


where v; is the stoichiometric coefficient of the reaction. 
Returning to reaction (Eq. 3.9) discussed before for which we have provided the 
thermodynamic data (Table 3.1): 


2H2(8) + O2(g) > 2H20(/) 


Like the entropy (A, S°) and enthalpy (A, H°) of reactions, the Gibbs free energy 
(A,G°) can be calculated from the product-minus-reactants (Eq. 4.6). 

A,G? = [2 x AG? ŒO) | z [2 x AG? (Ha(g)) + 1 x AG*(02(8))] = 
—474.26kJ, 

indicating that water is the stable phase at earth surface conditions. 

This Gibbs free energy change (A,.G°) can also be calculated based on the changes 
in entropy (A, S°) and enthalpy (A, H°) of reactions (Eq. 4.3): 

A,G° = A, H°—TA,S° = —571.66 x 1000—298.15 x —326.68 = —474.26 kJ 

and the results are identical. 

The first approach via the products-minus-reactants approach is simpler but is 
restricted to SATP conditions because AG‘ values are tabulated for SATP condi- 
tions, while the second approach via enthalpies and entropies can be used at other 
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temperatures as well. However, calculation of Gibbs free energy of reactions at non- 
standard temperatures via the second approach only holds if enthalpies and entropies 
do not vary significantly with temperature. In this introduction course, we make this 
assumption, being aware that both enthalpy and entropy vary with temperature. 

Changes in Gibbs free energy during chemical reactions can be used to perform 
work, to heat our homes or to support the metabolism of organisms. As an example 
of the latter, we calculate the A,G° for microbially mediated reactions. For aerobic 
(involving oxygen) and anaerobic (in the absence of oxygen, involving sulphate) 
degradation of model organic matter CHO (AG%. = —128.7 kJ mol~!) we have the 
following reactions: 


CH20(s) + O2(g) — CO2(aq) + H20(/) 


2CH20(s) + SOF (aq) > H2S(g) + 2HCO; (aq) 


The AG? values for these two reactions are — 494 and — 205 kJ, respectively, or 
normalized to one mole of CH20: — 494 and — 102 kJ. First, both aerobic respiration 
and sulphate reduction provide enough energy needed for growth and maintenance 
of the microbes involved (about — 10 kJ per mol). Second, about five times more 
energy is liberated during aerobic respiration than during sulphate reduction. Where 
is the missing energy? It is transferred to the reaction products: the hydrogen sulphide 
produced can react with oxygen: 


H2S(g) + 202(g) > H2S04(aq) 


and the A, G° value for this reaction is — 329 kJ. Microbial communities are highly 
efficient in using energy and the product of one organism is often the reactant for 
another organism, i.e., these reactions are coupled, with the result that most of the 
energy potentially available is eventually used. 


4.1 How Gibbs Free Energy Depends on Conditions 


The Gibbs free energy change during a process provides information whether a 
reaction or phase transformation will be spontaneous at constant temperature and 
pressure (AG < 0), but chemists and earth system scientists are often interested 
how the natural variables temperature or pressure impact the direction of processes 
or the stability of certain substances. We will therefore derive an equation, the funda- 
mental Gibbs equation, expressing the dependence of Gibbs energy on temperature 
and pressure. The term fundamental is used because P and T are the basic natural 
variables. We will use basic calculus for the derivations. 
Starting from the definition of the Gibbs free energy (Eq. 4.2), 
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G=H-TS 
and substituting the formal definition of H (Eq. 2.11): 
H=U+PYV, 
We arrive at, 
G=H-TS=U+PV-TS. 


Next, we take the total differential (i.e., the differential of terms of each variable 
in turn) 


dG =dU + PdV + VdP — TdS — SdT (4.7) 
Making use of the first law of thermodynamics: 
dU =q — PdV 


and the definition of entropy (Eq. 3.5): 


rewritten as TdS = q, 
we arrive at 


dU = TdS — PdV 
Substituting this equation into Eq. 4.7 leads to 
dG = TdS — PdV + PdV + VdP — TdS — SdT 
Simplifying by cancelling terms, we eventually arrive at the fundamental Gibbs 
equation which expresses the dependence of Gibbs free energy on changes in pressure 
and temperature: 


dG = VdP — SdT (4.8) 


Another way to quantify this dependency is to use the total differential (see Eq. 1. 
9 in Box 1): 


aG aG 
dG =(—) dp+(—) aT (4.9) 
IP); ðT)p 
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Comparing these equations, the partial derivative with respect to temperature is 


minus entropy 
( a) =-S (4.10) 
OT } p 


Hence, the entropy determines the temperature sensitivity of the Gibbs free 
energy. The larger the entropy, the greater the temperature dependence. Moreover, the 
Gibbs free energy of a compound declines (becomes more negative) with increasing 
temperature (because S values of compounds are always positive). 

Similarly, the partial derivative with respect to pressure is volume 


aG 
(5). =V (4.11) 


Consequently, the volume of a substance determines how its Gibbs free energy 
changes with pressure. Solids and liquids have small molar volumes (compared to 
gases) and the pressure dependence of the Gibbs free energy for solids and liquids 
is rather small. This is evidently not the case for gases. For an isothermal process, 
we can rewrite Eq. 4.8 as follows: 


nRT 
aG = var =( z Jap 


Integrating this equation to find AG 


Pr Pr (nRT Pry Py 
AG =f dG = / dP = nRT | —dP =nRTIn{ —) (4.12a) 
P; PA P p, P P; 


where we have used the ideal gas law to eliminate V and f A = Inx for the 
integration. 

This equation expresses the isothermal change in Gibbs free energy of a gas when 
pressures changes from P; to Pr. 

The change in Gibbs free energy (AG) at pressure P can then be expressed relative 
the standard Gibbs free energy change at | bar (AG”) 


P 
AG = AG? + nkTin( 2) = AG? +nRTIn(P) (4.12b) 


This latter equation is only valid if pressure is expressed in bar units, but the ratio 
(4) is dimensionless and Eq. 4.12a can be used when other units or other reference 


levels are adopted. 
Equation 4.12 has been derived for gases expressed in pressures, but an equivalent 
expression can be derived when considering ideal liquid or solid solutions. In this 
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case, the concentration [c] in mol L~! ([c] = y) is the analogue of pressure and the 
change in Gibbs free energy is: 


AG = AG? + nRTIn| £ | (13a) 
c 
where [c°] = 1 mol L~! for an ideal solution. For non-ideal solutions, we use 
thermodynamic activity a rather than concentration [c]: 


AG = AG? +nRTlna (13b) 


Thermodynamic activities for pure liquids and solids are | and the second term 
then becomes equal to zero and disappears as it should be because AG = AG? for 
pure phases. Later in the course (Chap. 5) we return to Eqs. 4.12 and 4.13 when we 
discuss equilibria of solutions. 


4.2 Phase Equilibria 


In the previous section we have derived the Gibbs free energy equation for a single 
phase and compound as a function of temperature and pressure. Now we are going 
to study the relationships between phases and how these depend on pressure and 
temperature conditions. A phase is a part of a system which is homogenous and 
separated from other phases by a definite boundary. As discussed earlier we deal 
with gas, liquid, and solid phases. Gases are always present as single phase because 
they mix well. A liquid containing a single component will only be one phase, but a 
mixture of liquid may be one phase when well mixed (ethanol in water) or multiphase 
when they do not mix (water—oil). Solids are often present in multiple phases, even 
a single component may be present in multiple phases (e.g., carbon in the form 
of graphite, diamond or fullerene, CaCO; as calcite or aragonite, Al,SiO5 as the 
minerals kyanite, sillimanite or andalusite). 
For phases A and B co-existing in equilibrium we can write: 


AG = 8G —dGg=0 (4.14) 


where dG 4 and dG are the Gibbs free energies of phases A and B. For both phases 
we can write the Gibbs free energy using the fundamental equation (Eq. 4.8) 


dGa = Vad P = SadT and dGpg = Vgd P — SpdT 
Since phases A and B are in equilibrium, dG 4 = dG g, hence 


VadP — SadT = VgdP — SgdT (4.15) 
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Which after re-arrangement reads 


VadP — VgdP = SadT — SgdT 


(Va — Vg)d P = (S4 — Sp)dT 


Introducing AV = (V4 — Vg) and AS= (S4 — Sg), we can rewrite this as 


AVdP = ASdT 

or 
dP AS aes 
dT AV ee 


Since we are at equilibrium AH = TAS we arrive at 


dP AH 
(4.16b) 


dT TAV 


This is the Clapeyron equation, and it describes the slope of the line between two 
phases in equilibrium (also called coexistence curve) in the pressure-temperature 
domain as a function of changes in molar volume (AV) and molar enthalpy (AH) 
or entropy (AS). These slopes are usually positive (Fig. 4.1) because entropy and 
volume changes are positive going from a solid to a liquid or liquid to a gas/vapour. 


Fig. 4.1 Generic one 
component phase diagram. 
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The Clapeyron equation can be used to determine if a phase transition is likely to 
take place under specific conditions. For instance, considering solids with a typical 
entropy of fusion (AS fusion) Of © 22 J mol! K~! and molar volume change during 
fusion (AV fusion) of © 4 x 10% m? mol! (i.e., 4 ml per mol), would yield 


dP AS 22 


ae AG ee 5.5 x 10*°PaK—! = 5Sbar K`! 
x 


Accordingly, a pressure increase of ~55 bar is required to change the melting 
temperature of solids by one degree. The Clapeyron equation is frequently used in 
solid-earth sciences because it relates the (geo)thermal gradient ( m = a) to changes 
in volume and enthalpy and provides then an estimate for the temperature T at which 
minerals are transformed or melting (solidus). 

This simple calculation can also be applied to water—ice equilibria although the 
volume change AV going from ice to liquid is negative (— 1.63 mL or — 1.63 x 1076 
m? for one mol of H20), indicating that the P to T slope is negative (Fig. 4.2). Another 
application is the calculation of freezing point increases with altitude (see Box 3). 

While volume changes for solid—liquid transition are rather limited (in the order 
of mL per mol), this is not the case for gases (about 24 L per mol) and slopes of 
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Fig. 4.2 Phase diagram of water. (Note the negative and steep slope of the ice-liquid water transition 
because of the negative volume change) 
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liquid—gas co-existence are flatter (Fig. 4.1). The change in volume of the system 
can thus be approximated with the volume of the gas, i.e. AV = V,,;. Assuming an 
ideal gas, we can rewrite Eq. 4.16 to 


dP AH AHP AHP 


dT TAV TRT RT? 


Next, we isolate pressure P on the left-hand side: 


dP  AHdT AHdT 
P RT RT 


(4.17) 


Making use of f a? = InP and integrating from P4, Tı to P2, T2 we arrive at the 
Clausius-Clapeyron equation: 


Pi AH(1 1 
n= = (4.18) 
P RAIL D 


The Clausius-Clapeyron equation is very useful for liquid—vapour/gas equilibria 
such as vapour pressure or boiling points calculation. As an example of the latter, 
we return to the Himalayas (Box 3). 


Box 3 Freezing and Boiling Points in the Himalayas 
Consider a mountain high in Himalaya with an atmospheric pressure of 35 kPa 
(altitude about 7.6 km) and calculate the freezing point increase. The molar 
enthalpy of melting (A H fusion) is 6.01 kJ mol7!, the volume change from ice 
to liquid is — 1.63 x 10~° m? for one mol of H>O and the pressure difference 
(AP) is 35 x10? — 1.013 x 10°= — 6.63 x 10* Pa. 

Using the Clapeyron equation (Eq. 4.16b), 


dP AH 6.01 x 103 


= = = —1.35 x 10’ PaK7! 
dT FAV 273.15-1.63 x 10-6 ae 


Hence: AT = ECT = 0.005 K. 

The freezing point is thus 0.005 °C higher. 

A similar calculation can be made for the lowering of the boiling point 
using the Clausius-Clapeyron equation (Eq. 4.18). The boiling point at 1 atm is 
100 °C but it will be lower high in the mountains because of lower atmospheric 
pressure. Using Eq. 4.18, with Pı = 1 atm =1.013 x 10° Pa and Tı = 100 °C 
= 273.15 K and the standard enthalpy of vaporization of water (A HĦyap) of 
40.7 kJ mol7!: 


4.3 Phase Diagrams 45 


1.013 x 10° 40,700 1 1 
n = — x 
35 x 10° 8.314 (373.15 Th 


we obtain a boiling point of 345 K (= 72 °C) at 7.6 km altitude. It takes much 
more time to boil an egg high in the mountains because the heat transfer from 
the water to your food is less. 


4.3 Phase Diagrams 


Materials exist in the solid phase at high pressure and/or low temperature conditions, 
but eventually transform into gas/ vapor at low pressures or high temperatures. Phase 
diagrams are a common way to depict which phase is present at a given P and T 
(Fig. 4.1). The phases (solid, liquid and gas) are separated by boundary lines, i.e., 
the lines represent conditions at which two phases co-exist. The Clapeyron relation 
(Eq. 4.16): af = aS gives the slope of these lines. 

There are three types of lines: (1) the line between the solid and liquid phase 
shows how the melting points (temperature of solid to liquid transformation) change 
with pressure, (2) the line between liquid and vapor phase shows the dependence 
of boiling points (temperature of liquid to vapour transformation) on pressure, and 
(3) the line between solid and vapor delineates the transition from a solid to gas 
(sublimation and condensation). The triple point represents the P-T combinations 
at which the three phases co-exist in equilibrium and the critical point represents 
a combination of temperature and pressure beyond which a gas cannot be liquified 
regardless of the pressures and temperature. The state of matter beyond the critical 
point is neither a gas or liquid; it is a supercritical fluid that possesses some typical 
liquid properties such as high density and ability to act as solvent, but the viscosity 
and diffusion coefficients are more alike in a gas. 

The ‘normal’ melting point is where a horizontal line at 1 bar or 1 atm crosses the 
solid—liquid equilibrium lines, and the ‘normal’ boiling point is crossed at the liquid— 
gas (vapor) equilibrium lines (Fig. 4.2). These are called normal melting/boiling 
points because they correspond to Earth-surface conditions, e.g., 0 and 100 °C for 
the normal melting and boiling points for water. In Box 3/Sect. 4.2 we have calculated 
how boiling points drop when pressure drops or increases when pressure increases. 

The phase boundary lines constrain the freedom of selecting pressure and temper- 
ature conditions for a system to be in thermodynamic equilibrium. Within a single 
phase, one can freely vary P and T independently. At the boundary between two 
phases, i.e. at the equilibrium lines, either P or T can be chosen, the other one is 
then fixed. Gibbs has formalized the number of independent variables (= degrees of 
freedom) to describe a multicomponent system. Gibbs’ phase rule states that 


F=C-P+2 (4.19) 
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where F is degrees of freedom, C is the number of components and P is the number 
of phases in equilibrium. The + 2 term relates to temperature and pressure. Consider 
a phase diagram for one single component, e.g. H2O (C=1) (Fig. 4.2). Within fluid 
water, or within the domain of ice or water vapor, pressure and temperature can be 
combined in multiple ways because there are two degrees of freedom (F = 1-1 + 
2 = 2). At the interface of water—ice, there is only one degree of freedom (F = 1-2 
+ 2 = 1), i.e., if temperature is chosen, pressure is fixed or the other way around. 
The boundary lines between mineral phases are therefore called univariant curves in 
petrology and mineralogy. At the triple point, ice, water and water vapour co-exist 
(water is freezing and boiling at the same time) and Gibb’s phase rule indicates that 
there are no degrees of freedom (F = 1-3 + 2 = 0). In other words, the triple point 
has a fixed pressure and temperature for a component, e.g., for water it is 0.0098 °C 
and 6 x 10° atm. 

Phase diagrams are instructive to show that melting (solid to liquid) can occur 
either through increasing the temperature or through lowering the pressure (except 
for water because of its negative slope). Similarly, boiling of a liquid can not only 
occur by increasing the temperature, but also induced by lowering the pressure, i.e., 
it is possible to boil cold water under low pressure conditions. From the quantitative 
treatment in Sect. 4.2, it is clear that solid—liquid lines are normally steep because 
the volume changes involved are small, while liquid—gas equilibrium lines are rather 
flat because volume changes are large (about three orders of magnitude for ideal 
gases). Moreover, the liquid—gas equilibrium lines usually show substantial curva- 
ture because of the temperature dependence of enthalpies of vaporization. Finally, 
this course is limited to single component phase diagram, but the same principles, 
including Gibbs’ phase rule, apply to multicomponent systems. However, it can get 
rather complex. 
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Part II 
Equilibria: Solutions, Minerals, Acid-Base 
and Redox Reactions 


Chapter 5 A) 
Introduction to Equilibrium get 


Abstract This chapter defines equilibrium in mixtures and ideal solutions and intro- 
duces the equilibrium constant, its relationship with the Gibbs free energy and its 
dependence on temperature. Homogenous and heterogenous equilibria are distin- 
guished and solid—gas, liquid-gas and mineral-solution equilibria are presented, 
including Henry’s law and the solubility product. 


Keywords Equilibrium constant - Reaction quotient - Solubility product - van ‘t 
Hoff equation 


Our thermodynamic classes have provided us the tools to determine whether a reac- 
tion at Earth surface conditions is spontaneous or not, i.e., whether it has a negative or 
positive Gibbs free energy. For instance, the dissolution of sodium chloride (kitchen 
salt) is an endothermic process (consuming heat and cooling the water; A; H° > 0) 
but occurs spontaneously because it has a negative Gibbs free energy (A, G°) of —9 
kJ mol~!. Spontaneity of a reaction does not imply that the reaction goes to comple- 
tion, i.e., equilibrium is not where all reactants are converted into products. While 
addition of kitchen salt to water initially indeed results in the dissolution of NaCl, no 
further dissolution is observed after that about 357 gr has been added to | kg of water. 
If more than 357 gr of NaCl were dissolved in water, it would precipitate (a backward 
reaction). The system water-NaCl has reached equilibrium, i.e., the concentrations 
of reactants and products remain constant over time. An equilibrium state is attained 
as the rate of forward (from reactants to products) and backward (from products to 
reactants) are equal and there are no further changes in concentrations. In System 
Earth many reactions can be fully described, or be approximated, by an equilibrium 
approach: the exchange of gases between water and air, proton transfers (acid-base 
chemistry) and electron transfers (redox reactions) and the precipitation and disso- 
lution of minerals. For instance, the dissolution of calcium carbonate leads to karst 
processes and the precipitation of calcium carbonate causes formation of stalactites 
and stalagmites in caves. 

Itis therefore useful to deepen our understanding and further develop our quantita- 
tive tools for predicting equilibria. In this introduction, we will focus on ideal solution 
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chemistry and equilibria at Earth-surface conditions and leave non-ideal solutions 
and high-temperature and high-pressure equilibria during formation of magmatic and 
metamorphic rocks to subsequent courses. Following the introduction of the equi- 
librium constant and its relation to the Gibbs free energy, we will focus on equilibria 
in ideal solutions involving minerals, acids and bases, and redox reactions. 


5.1 The Equilibrium Constant and Its Relation 
to the Gibbs Free Energy 


Nitrogen dioxide (NO2) is an atmospheric gas resulting from combustion engines 
(cars, buses, trucks). It is a deep brown gas that contributes to smog above cities. In 
the atmosphere, it is in equilibrium with dinitrogen tetroxide (a colorless gas): 


2NO2(g) <> N204(8) (5.1) 


The transformation of NO2(g) to N2O4(g) causes the color to fade until an equilib- 
rium is established, and this reaction has therefore been well studied in the laboratory. 
Irrespective of the starting conditions (only NO2(g), only N2O4(g), high NO2(g) with 
low N20; (g), or the other way around), the final composition of the reaction mixture 
is always the same when expressed as the ratio of product to reactant concentrations 
(the law of mass action). Specifically, in this case the calculated concentration ratio 
moo j is about 216 at 25 °C and 1 bar, irrespective of starting conditions. 

To generalize this phenomenon, we introduce the reaction of œ moles A and 6 
moles B to form y moles C and ô moles D 


aA+BB>yC+6D (5.2) 
and define a reaction quotient Q: 


Cc’ D° 
~ A® BB 


(5.3) 


This reaction quotient is dimensionless and general, i.e., it can always be calcu- 
lated irrespective whether there is an equilibrium or not. When the system is at 
equilibrium, the value of Q is called the equilibrium constant K: 

Y pê 
= Ceq D eq 
~ B 

Ae, Bea 


(5.4) 


The value of K expresses the composition of the system, i.e., concentrations 
or partial pressures of the reactants and products at equilibrium (taking reaction 
stoichiometry into account) and has a constant value at a certain temperature and 
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pressure. The equilibrium constant has no unit because they cancel out. Reversing a 
reaction causes an inversion of K, and multiplying the coefficients (a, B, y, 8) by a 
common factor raises the equilibrium constant to the corresponding factor. If the K 
value is very high, the equilibrium mixture will consist mainly of products C and D. If 
K < 1, then reactants will dominate the mixture, while similar quantities of reactant 
and products are expected if K is around 1, depending on the reaction stoichiometry 
of course. 


K >> 1 Products dominant 
K %1 Similar quantities 
K << 1 Reactants dominant 


Comparing the reaction quotient Q and the equilibrium constant K provides infor- 
mation on the direction of change, i.e., whether a process is spontaneous. If Q > K, 
then the net reaction goes from right to left (products to reactant, backward reaction 
is spontaneous) and if Q < K then the reaction goes from left to right (reactants are 
transformed into products and the forward reaction is spontaneous). If Q = K then 
the reaction is in equilibrium. 


Q > K backward reaction is spontaneous 
Q < K forward reaction is spontaneous 
Q = K reaction is in equilibrium 


In thermodynamics lectures we have seen that spontaneity is related to the Gibbs 
free energy change at constant pressure and temperature; there should thus be a 
relation between the reaction quotient Q, the equilibrium constant K, and the sign of 
the Gibbs free energy change (A,G). 

Returning to reaction (5.1), the Gibbs free energy change for the gas reaction 
(A,-G) of NO2(g) to N204(2) 


A-G = AG(N204(g)) — 2 x AG(NO3(g)) (5.5) 


We use A,G rather than A, G° because the reaction is not occurring under standard 
conditions (i.e., for a pure gas). During the thermodynamic classes (Eq. 4.12) we have 
derived that 


P 
AG = AG? + arin( $) (5.6) 
0 


where R is the universal gas constant, T is absolute temperature, AG” is the standard 
Gibbs free energy change for a pure gas, P, is the partial pressure of gas g, and Po is 
the total pressure of | bar, respectively. Combining Eqs. (5.5 and 5.6), we then arrive 
at 
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Using A,G? = AG* (N204(g)) — 2 x AG% 4 (NO2(g)), we can rewrite Eq. (5.7) 
as follows, 


Pxyn04 
P (78) 
A.G = A,G° + RT Inf 2%) — 2RT in( X82) = A,G° + RT In 
Po Po (2) 
Po 


(5.8) 


Accordingly, the Gibbs free energy of reaction A, G contains the standard Gibbs 
free energy change of reaction (A, A that is independent of the partial pressures 


ae: 


of the two gases and a term RT In 5 that depends on the partial pressure of the 


reactant and product. 


(Ce) 


nov? ÍS the reaction quotient Q for reaction 5.1, if expressed 
Po 
in partial pressures and a total pressure Po of | bar. 


Generalizing this to any reaction, we obtain 


Moreover, this ratio 


— 


A,G = A,G? + RT InQ (5.9) 


stating that the Gibbs free energy change of a reaction is the sum of the standard 
Gibbs free energy change of the reaction (A, G°), calculated from the Gibbs free 
energy of formations, and the term RT/nQ, where the latter provides a correction 
term for the actual composition of the mixture. 

As an example, consider reaction 5.1 and a reaction vessel at 298 K with 0.350 
bar NOx) and 0.65 bar N2O4(g), The standard Gibbs free energy of formations of 
the two gases are AG" 51.3 and 99.8 kJ mol“, respectively. 


A,G=A,G? +RTO 


0.65 
= 99,800 — 2 x 51,300 + 8.314 x 298 In 


5 ast = 13303 


The A,G 4 0, the system is thus not at equilibrium. Specifically, the net reaction 
proceeds towards reactants, i.e., NO2(g), because A,G > 0. (The same conclusion 
could be reached using the reaction quotient Q and the equilibrium constant K). 

At equilibrium, the reaction quotient is equal to the equilibrium constant (Q 
= K), and the Gibbs free energy of the reaction should be zero (A;G = 0; see 
thermodynamic lectures, Eq. (4.4)). Consequently, 
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0=A,G=A,G° + RTInK 


Hence, 
A,G° = —RT In K (5.10a) 
or alternatively written, 
K = exo( =) (5.10b) 
RT 


Equations (5.10a and 5.10b) are completely general and refer to all types of 
reactions. They provide a direct link between thermodynamic data and theory on 
the one hand and the composition of equilibrium mixtures on the other. However, 
Eqs. (5.10a and 5.10b) implies that the mixture composition is expressed in activities 
rather than concentrations (see Eq. 4.13). In this course, we restrict ourselves to ideal 
solutions, i.e., ion pairing and interactions between ions and the solvent, are ignored 
and we express the equilibrium constants and reaction quotients in concentrations 
[c], rather than thermodynamic activities a, where [c] = Z for a gas (Po = 1 bar), 


[c]= Sa for a solute (Cp = 1 mol L7!) and [c] = 1 for solids or the liquid medium. 
Activities (or effective concentrations) and concentrations are related via activity 
coefficients. However, the use of concentrations rather than activities implies the 
need of units for equilibrium constants. 

Equilibrium mixture compositions are temperature dependent and there are two 
ways to correct equilibrium constants for lower or higher temperatures. The first 
method uses Eq. (5.10a and 5.10b) to link the equilibrium constant K with the A, G° at 
the temperature of interest calculated from the entropy (A^, S°) and enthalpy (A, H°) 
of reactions (Eq. 4.3 in thermodynamics). Alternatively, the van ‘t Hoff equation is 
used. 

Combining, Eq. (5.10a): A-G? = —RT In K. 

and Eq. (4.3): A,G° = A,H° — TA,S°. 

we obtain the van ‘t Hoff equation after rearrangement, 


A,S° A, H® 


In K = (5.11) 
R RT 
Or in its more useful form going from T; to T3, 
In Kọ = In K pa ar (5.12) 
on Ne ae 


where Tı is usually 298 K (SATP). If the reaction is exothermic (A,;H° < 0), 
then K decreases with increasing temperature and equilibrium shifts towards more 
reactants, i.e., to the left. Conversely, equilibria shift to the right and K increases with 
temperature for an endothermic reaction (A, H° > 0). 
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5.2 Heterogenous Equilibria, Including Mineral Solubility 


Equilibrium can be established when all reactants and products are in one single phase 
(homogenous equilibrium), or multiple phases are involved (heterogenous equilib- 
rium). Acid-base reactions in solution are an example of the former, while mineral 
dissolution/precipitation reactions involving a solid and a liquid phase are heteroge- 
nous. This section will present three types of heterogenous equilibria (solid—gas, 
liquid—gas and mineral-solution). Homogenous acid-base equilibria will be covered 
in the next chapter. 

Solid—gas equilibria: The production of lime (CaO, s) from limestone (CaCO3, s) 


CaCO3(s) <> CaO(s) + CO2(g) 


is a heterogenous reaction with an equilibrium constant K defined as 


k= [COICO _ cg, 
[CaCO3] 


Note that [...] are used to indicate concentrations rather than activities and that 
solid phase and liquid can be excluded from heterogenous equilibria mass action laws 
since they have an (activity) value of one. Consequently, the equilibrium concentra- 
tion of carbon dioxide gas in a system with both solid phase CaO and CaCO; is 
independent of the amount of these solid phases, as long as both are present. 

Liquid—gas equilibria: The solubility of gases in water is another example of 
heterogenous equilibrium. The solubility of a gas in water is described using Henry’s 
law: 


Cgas = Ky: Poas (5.13) 


where Ceas is the equilibrium concentration of a gas in water (mol kg), Pgas is 
the partial pressure (atm) and Ky is the Henry’s law constant (mol kg~! atm7'), 
Henry’s law constant is a function of temperature and salinity of the water. Consider 
the solubility of oxygen in water: 


O2(g) + H2001) <+> O2(aq) 


The corresponding equilibrium constant is Kyo, = new, 


Figure 5.1 shows the concentrations of dissolved oxygen as a function of temper- 
ature for freshwater and ocean water. Oxygen and most other gases dissolve better in 
freshwater than seawater because of a salting out effect. Dissolved oxygen concen- 
trations in water decline with warming of the water. Cold waters in equilibrium with 
the atmosphere therefore have higher oxygen contents than warm waters. 
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Fig. 5.1 Solubility of oxygen as a function of temperature in freshwater and seawater 


Mineral-solution equilibria: The availability of certain minerals constrains the 
composition of natural waters and vice versa the composition of solutions can deter- 
mine whether minerals dissolve or precipitate. Mineral equilibria are written as disso- 
lution reactions (i.e., mineral on the left-hand side and ions on the right-hand side) and 
characterized by a solubility product. The solubility product of the mineral fluorite 
(CaF2,,)) is based on the following reaction: 


CaF2(s) <> Cat (aq) + 2F (aq) (5.14) 


with the following equilibrium constant 


[c+] [F]; 
[CaF>] 


K = 


where the subscript eq indicates that these are equilibrium concentrations. Recalling 
that solid phases have a (molar fraction) value of one, we can then define the solubility 
product for fluorite as 


2 
eq 


Ksp = [Cat] [F] (5.15) 


The solubility product K,, is the equilibrium constant (K) for the special case 
that a solid phase (mineral) is dissolving. This can be generalized to 


Cr Ay <> xC + yAee 
and 


Ka = [L |aE (5.16) 
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Returning to the solubility of fluorite (reaction 5.14), if we know the equilibrium 
concentrations of dissolved calcium (2 x 10“ M) and fluoride (1.6 x 10+ M), we 
can use Eq. (5.15) to calculate the solubility product as 5.1 x 107!?. A solubility 
product is often reported as a pK, value, i.e. the —log9(K sy), in this case pK sp = 
11.29. 

The same result would be obtained if calculated using Eqs. (5.10a and 5.10b) and 
thermochemical data with AG* values of —1175.6, —553.6, —278.8 kJ mol`! for 
fluorite, dissolved calcium, and fluoride, respectively. First, one calculates the Gibbs 
free energy of reaction (A,G°) :Type equation here. 


A,G? =[1x A G% (Ca**(aq)) + 2 x A G4 (F (aq) ] 
— [A G4 (CaFy(s))] = 64.4 kJ 


and then uses Eqs. (5.10a and 5.10b): 


K ao = 5x 10°" K 11.29 
= ex = ex =5x > =11. 
PLR P\ 8:31. 298 cee 


Accordingly, solubility products link thermochemistry with solution chemistry 
for solutions in equilibrium. 

Solubility products can also be used to constrain the composition of natural waters. 
Consider a pool in equilibrium with the minerals fluorite (CaF2, pK,, = 11.29), 
calcite (CaCO3, pK,, = 8.48) and dolomite (CaMg(CO3)2, pK, = 19), knowing 
that the dissolved fluoride concentration is 10% M. From the fluorite equilibrium, 
we then obtain a dissolved calcium concentration of 5.13 x 10% M, next the calcite 
equilibrium can be used to derive the carbonate ion concentration 6.46 x 10 M 
and finally the dissolved magnesium concentration is calculated making use of the 
solubility product of dolomite (4.68 x 10 M). This example clearly shows that 
sharing common ions has major implications for the solubility of minerals. 

The precipitation and dissolution of minerals depends on the saturation state of 
that mineral. For the dissolution of calcium carbonate in the form of calcite: 


CaCO3(s) <> Ca?* (aq) + COF (aq) (5.17) 


we first define the ion product (IP), which is a special form of the reaction quotient 
Q for mineral equilibria. 


IP = [Ca”*][Co3" (5.18) 


with, in this case, a unit of mol? kg’, and then the degree of saturation or saturation 
state (Q2) 


[eico] _ [oxt*J[cor] _ rp 
5 fear], [COF], ~ 7 o 


q eq 
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Calcite will dissolve if Q < 1, i.e., in undersaturated water, and precipitate from 
solution if Q > 1, i.e., if waters are supersaturated with respect to the mineral. 
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Chapter 6 A) 
Acid-Base Equilibria get 


Abstract This chapter presents acid-base equilibria and simple methods to calculate 
the pH of solutions. The inorganic carbon system in water is introduced, including 
the concept of alkalinity, and methods to solve carbon dioxide equilibria in water 
are discussed. Earth system science relevant examples such as rainwater, surface 
waters in equilibrium with calcium carbonate minerals, soda lakes and the physical 
chemistry of karst systems are presented. 


Keywords Acid - Base - Dissolved inorganic carbon - Alkalinity - Karst © Soda 
lake + Rainwater - Ocean acidification 


6.1 Introduction 


Although protons are produced and consumed in many chemical processes, most 
natural waters fall within a narrow range of pH values of 6-9 because of interactions 
with minerals and transfer of protons (i.e., hydrogen ions) among acids and bases 
in solution. Most proton transfer reactions in solution are very fast (i.e., reach an 
equilibrium in milliseconds) and a (homogenous) equilibrium description is thus 
appropriate. A proton does not exist as such in water because it reacts fast and 
strongly with water to form a hydrated proton (hydronium), H3O0*, which in turn 
associates with hydrogen bonds to additional water molecules. It is macroscopically 
not possible to distinguish between the various hydrated proton species and we will 
therefore use H* and H30* interchangeably and refer to H* as the proton or hydrogen 
ion in solution. 

An acid is a substance that dissociates in water releasing protons into solution, 
i.e., it is a proton donor, while a base accepts protons. Chemical species whose 
formulas differ only by one hydrogen ion are called conjugate acid-base pairs. An 
acid HA dissolved in water dissociates into a proton and a base A“: 


HA(aq) + H20(1) <> H30* (aq) + A~ (aq) (6.1) 
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for which we can write the general equilibrium constant K 


[Hs0*][A7] 
[HA][H20] 


Noting that the liquid medium has a (activity) value of 1, we can then define the 
acid dissociation or acidity constant 


Ht |[A7 
a= uae (6.2) 


A similar treatment can be followed for base B added to water 
B(aq) + H,0(1) + OH” (aq) + BH? (aq) (6.3) 
with the base dissociation or basicity constant 


[oH |[BH*] 


K, = iB] ; 


(6.4) 


but we will elaborate it from the acid side of the story. 

The extent of dissociation, i.e., the proton release by the acid HA and its transfer 
to water acting as the base in Reaction 6.1, depends on the strength of the acid, i.e., 
the value of K,. Strong acids nearly fully dissociate into H* and A~ and almost no 
HA remains in solution. Weak acids dissociated partially and after equilibration HA, 
H* and A” are found in solution. Acids stronger than H30* are called strong acids 
and their conjugated bases are very weak. Weak acids are less strong than H3O* but 
stronger than H20 (Table 6.1) and their conjugated bases are also weak. 

Reactions 6.1 and 6.3 show that water can act as a proton donor as well as proton 
acceptor. The dissociation or self-ionization of water, 


2H20(1) <> H30+ (aq) + OH” (aq) (6.5) 


leads to the ion-product constant for water, i.e., the equilibrium constant for self- 
ionisation of water, 


K = [H30*][OH ] = 10-* (6.6) 


Some acids contain multiple protons, and these are called polyprotic acids (e.g. 
carbonic acid, phosphoric acid). These polyprotic acids dissociate stepwise, and each 
dissociation step is characterized by its own acid dissociation constant. In solutions 
containing multiple acids and bases, protons are transferred from the stronger acid 
to the stronger base to yield the weaker acid and base. 

For notational simplicity, we will drop the phase descriptions aq for solutes and / 
for liquid water, but the s for solids is retained. Acid and base dissociation constants 
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Table 6.1 Acid and base dissociation constants of conjugated acid-base pairs 
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Name Acid pKa Base pKb 
Hydrochloric acid HCl Strong Cl” Very weak 
Sulfuric acid H2SO4 Strong HSO% Very weak 
Nitric acid HNO3 Strong NO37 Very weak 
Proton H30* H20 

Chromic acid H2CrO4 0.74 HCrO4 13.26 
Iodic acid HIO3 0.78 1037 13.22 
Sulfurous acid H2SO3 1.85 HSO37 12.15 
Bisulfate HSO47 1.99 S04? 12.01 
Phosphoric acid H3PO4 2.16 H2 P0O4? 11.84 
Hydrofluoric acid HF 3:2 F- 10.8 
Nitrous acid HNO2 3.25 NO27 10.75 
Hydrogen selenide H2Se 3.89 HSe~ 10.11 
Acetic acid CH3CO2H 4.76 CH3CO2 9.24 
Carbonic acid H2CO3 6.35 HCO37 7.65 
Bichromate HCrO47 6.49 CrO,? 7.51 
Hydrogen sulfide H2S 7.05 HS~ 6.95 
Bisulfite HSO3— 7.2 SO37- 6.8 
Dihydrogenphosphate H2PO4?- 7.21 HPO,?- 6.79 
Boric acid H3BO3 9.27 H2BO3 4.73 
Ammonium NH,t 9.3 NH3 47 
Bicarbonate HCO3~ 10.33 CO37— 3.67 
Biselenide HSe~ 11 Se?~ 3 
Biphosphate HPO, 12.32 PO, 1.68 
Water H20 OH 

Ammonia NH3 Very weak NH2~ Strong 
Hydroxy ion OH” Very weak o% Strong 


are often presented as pK values (— log K), and [H+] and [OH | concentrations 
are traditionally presented as pH and pOH values (— log [H+] and — log [OH™]). 
The pH scale is logarithmic implying that a 0.3 increase corresponds to a halving of 
concentrations (— log(2) = 0.3) and that the absolute change depends on the value: 
from pH 5 to 5.3 corresponds to a 5 yM proton concentration decrease from 10 to 5 
uM, while a pH change from 6 to 6.3 relates to 0.5 uM change in [H+] from 1 to 


0.5 uM. 


pH calculations: When calculating the pH of a solution in which an acid is 
dissolving, the first step is to evaluate whether the acid is strong (pKa < 0) or weak. 
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A strong acid will fully dissociate, i.e., all the acid HA added to the solution will 
be converted into H* and A~ and the eventual pH = — log [HA]. A 0.1 M solution 
with the strong acid HC] will thus have a pH of 1. 

The pH of a solution of a weak acid HA (with concentration C,) dissolving in 
pure water can be calculated with a (R)ICE table. 


(Reaction) [HA] [a+] [a7] 
Initial Ca ~o a 
Change =x x 
Equilibrium Ca-x X = 


The first row is for the reaction: HA <> H* + A7 

The second row specifies the initial conditions: the concentration of the acid added 
(Ca) that is initially undissociated so that [a7] is zero and [H+] is close to zero (1077 
for pure water). 

The third row gives the change due to reaction; x amount of [HA] is consumed 
and x [H+] and x [a7] are produced. 

The fourth row, equilibration, is simply the sum of rows two and three. 

Next, the acid dissociation constant (K4) and the concentrations in the equilibrium 
line in the (R)ICE table are combined in Eq. 6.2. 


_ [H*][A7] O Xx 


[HA] Cy—x 


which can be rewritten as a quadratic equation: 
x? + Kax — K,Ca = 0 


which can be solved for x to obtain the proton concentration, and thus the pH = — 
log(x). If Ca > 1000 times Ka, then x is very small relative to Ca, and Ca — x ~ Ca. 
The proton concentration can then be calculated from [H+] = VK,Ca. 

As an example, calculate the pH of 0.1 M Acetic acid solution with a pK, value 
of 4.76. The corresponding (R)ICE table reads 


(Reaction) [HAc] [Ht] [Ac] 
Initial 0.1 ~0 0 
Change =% X x 
Equilibrium 0.1 — x x x 


and the resulting pH = 2.88. 
Common ions, i.e., substances that are both in the acid/base system as well as 
in the solution, shift the equilibrium and thus the pH of the solution. Consider the 
above example, but now with 0.05 M Na-acetate as well in the solution. Na-acetate 
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has a very high dissociation constant and is fully dissociated, i.e., it is present in 
the form of acetate [Ac] and sodium ion [N at]. Sodium ions are not involved in 
proton transfer reaction and can be ignored. The (R)ICE table must be modified for 
the initial conditions to accommodate for the presence of acetate ions. 


(Reaction) [HAc] [Ht] [Ac] 
Initial 0.1 ~0 0.05 
Change =% x x 
Equilibrium 0.1 —x x 0.05 + x 


The corresponding equilibrium 


[H+] [Ac] x- (œ + 0.05) 


[HA] ` O0l-x 


and quadratic equation 
x? + (0.05 + Ka)x — 0.1 Ka = 0 


results in a pH of 4.46. The addition of the common ion acetate shifts the equilibrium 
to the left and the proton concentration declines (and pH increases). 

Buffer solutions: Solutions containing weak acids with their conjugated bases are 
buffer solutions because they resist drastic changes in pH. Quantifying the buffering 
capacity of natural waters containing multiple acid-base systems can be quite cumber- 
some. Therefore, we restrict ourselves to the Henderson-Hasselbalch equation for a 
single acid-base system to illustrate the principle. 

Starting with re-arranging Eq. 6.2 


WA 
[HA] 


to isolate the proton concentration on the left-hand side 


1 A] tA] 


[H+]  Ka[HA] Ka [HA] 


followed by taking the logarithms 


[A] _ [A] 
maj Re Oe A (6.7) 


pH = — log K, + log 

we have derived the Henderson-Hasselbalch equation. 
Dissolution of carbon dioxide in water results in the formation of carbonic acid 
and this provides most of the buffering of natural waters. In the next section we will 
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elaborate the carbon dioxide-water system in quite some detail, here we illustrate the 
buffering provided by the carbonic acid-bicarbonate system: 


H,CO; <> Ht + HCO, 
Consider a system with initial HCO; and HCO3~ concentrations of 1 mM each 


and a pH of 6.35, i.e., the pK, value. Using the Henderson-Hasselbalch equation, we 
obtain 


H = pK, +1 [meo] 6.35 +1 Ler] 6.35 
= a og ——— = 0. (0) = 0. 
aaa  [H2CO;] STOA] 


Addition of 10 u mol strong acid to 1 liter would lower HCO3~ and increase 
H2C0O; concentrations, but hardly impact the pH 


[10-3 — 1075] 


e—a _ 6.34 
[10-3 + 10-5] ae 


pH = 6.35 + log 


The ApH is 0.01 at the initial pH of 6.35, but it would be ~ 0.03 at a pH of 5.35 
and © 0.31 at a pH of 4.35. Accordingly, a buffer is optimal close to its pK, value, but 
still limits proton concentration changes to within a factor of two at pH values two 
units higher or lower than the pK,. Carbonic acid is a diprotic acid, i.e. a polyprotic 
acid with two protons, with pK, values of 6.35 and 10.33. It thus has the potential to 
buffer natural waters over a wide range of pH values. 
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Carbon dioxide is a gas that exchanges between water and air. When dissolved in 
water, it reacts with water to form carbonic acid, which in turn dissociates into 
bicarbonate and carbonate ions. This carbon-dioxide-water system provides much 
of the buffering capacity for natural waters, governs the pH of most natural waters, is 
pivotal to precipitation and dissolution of carbonate minerals, and governs the uptake 
of (anthropogenic) carbon by the ocean. 

The carbon-dioxide-water system is described by a few equilibria. Carbon dioxide 
gas dissolves in water and will then react with water to form carbonic acid (H2CO3) 
according to 


CO2(aq) + H20 < H2CO3(aq) (6.8) 
This equilibrium lies rather far to the left, hence most of the dissolved carbon 


dioxide remains in the form of CO2(aq), and we can analytically not distinguish 
between H2CO3(aq) and CO2(aq); these two pools are therefore lumped together 
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and termed carbonic acid (H2CO3) in this course. Other treatments of carbon-dioxide- 
water system sometimes use the terms HCO} or CO3 for the lumped pool. Accord- 
ingly, we can re-write Eq. 6.8 to directly link atmospheric carbon dioxide and carbonic 
acid as follows: 


CO2(¢) + H20 < HCO; (6.9) 
with the equilibrium expression based on Henry’s law (Eq. 5.13): 

[H2CO3] = Ky - Pco, (6.10) 
where Pco, is the partial pressure of carbon dioxide in equilibrium with water and 
Ky is Henry’s Law constant (mol kg~! atm~!) with a pKy value of 1.47 at 25 °C 
in pure water. Carbonic acid is a weak acid that partly dissociates into a proton and 
bicarbonate (HCO, ), the latter being another weak acid that dissociates into another 
proton and carbonate (C037). The relevant reactions are 

H2C0; <+ H* + HCO; (6.11) 
and 


HCO; < H* + C037 (6.12) 


with the corresponding equilibria and equilibrium constants (in freshwater and at 
25 °C) 


_ [EĦ][BCO;] as 


[MCO] oe 
K> = [H*][co;] = 1071033 (6 14) 
[HCO; | l 


Accordingly, the carbon-dioxide-water system is characterized by five aqueous 
species: [H+], [OH], [H2CO3], [HCO; ] and [Coz] that are linked via three 
equilibrium relations (the self- ionization of water, Eq. 6.6, and the first and second 
acid-base equilibria, Eqs. 6.13 and 6.14) with known constants. To calculate the 
abundance of the five species, we must solve the system and for this we need more 
information, specifically, two more relations. 

Various types of information can be used. These can either be information on 
the measured or known concentration of the aqueous species (e.g., known pH), 
or additional relations. For a system open to exchange with the air, Henrys’ law 
(Eq. 6.10) provides a link of carbonic acid (H2CO3) to the partial pressure of carbon 
dioxide in the air (Pco,). Solubility equilibria involving carbonate minerals can 
constrain the carbonate ion concentration. Natural waters are uncharged, and we can 
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thus use the charge balance equation as an additional constraint: the positive charge 
of protons should be balanced by the negative charge of hydroxide, bicarbonate and 
carbonate ions. 


[H*] = [OH] + [HCO;] + 2[CO;] (6.15) 


Note that the carbonate ion is counted twice in a charge balance because of 
its double charge. Alternatively, one can define a proton balance equation, a mass 
balance for protons 


[H*] = ES ae + [H+] inco, (6.16a) 
or its equivalent 
[H*] = [0H] + [HCO; ] + 2[c0; ], (6.16b) 


This proton conservation equation balances excess protons on the left-hand side 
with the recipe on the right-hand side and is in this case identical to the charge balance 
(Eq. 6.15). However, the proton conservation and charge balance equations can be 
different in more complex solutions such as seawater. 

Unfortunately, not all species of the system can be measured directly, but there is 
also no need because they are interlinked via equilibrium, mass balance and charge 
conservation equations. The four parameters most often measured are: the pH, the 
partial pressure of (Pco,), the dissolved inorganic carbon content and the titration 
alkalinity. The total dissolved inorganic carbon content, abbreviated as Cr or DIC 
(dissolved inorganic carbon) is defined as: 


DIC = Cr = [H2CO3] + [HCO; ] + [COZ]. (6.17) 


The titration alkalinity is the excess of proton acceptors over donors of a solution 
and is normally derived from an acidimetric titration. Alkalinity (TA) for our simple 
system is defined as: 


TA = [OH | +[HCO;] + 2[CO;] - [B+] (6.18) 


This alkalinity definition is rooted in the charge and proton balances presented 
earlier (Eqs. 6.15 and 6.16). Carbonate alkalinity (CA = [HCO;] + 2 [CO3 
predominates alkalinity of natural waters because of the presence of carbon dioxide 
in the atmosphere. 
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The previous section has provided the tools to solve carbon dioxide equilibria. 
However, solving carbon dioxide equilibria in natural waters, particularly in seawater, 
can be challenging because of the presence of multiple acid-base systems and 
other species interacting with the carbonate system. Fortunately, there are multiple 
computer programs available to assist us. Below we present a few simple, Earth- 
sciences relevant, cases that are instructive and can be done without such computer 
programs. First, it is important to identify whether the system is open (to an atmo- 
sphere with known composition) or closed (when interactions with the atmosphere 
are negligible). The open system approach applies to surface waters of streams, rivers, 
lakes, reservoirs, estuaries, and the surface ocean, while a closed approach is more 
appropriate for groundwater and subsurface waters in lakes and in the ocean inte- 
rior that are isolated from the atmosphere. Second, calculations for problems with 
known pH are usually simpler than those without a priori knowledge on pH, and we 
therefore start with the former. 

Case I The pH is known 

This type of problems is most conveniently solved using the ionization fraction 
approach. The derivation of the ionization fraction starts with rewriting the first 
carbonic acid equilibrium (Eq. 6.13) to isolate bicarbonate on the left-hand side: 


[E Reo] -4 _ Kı[H2C0;] 
Ki = co) [HCO; ] = m (6.19) 


Next, we substitute this equation in the second carbonic acid equilibrium (Eq. 6.14) 


(6.20) 


and we obtain after re-arrangement an expression with the carbonate ion isolated on 
the left-hand side: 


__ KıK[HCO;] 


(A M 
[co] ] = mT (6.21) 


Combining these equations with the definition of the dissolved inorganic carbon 
content (Eq. 6.17) 


DIC = [H2C0;] + [HCO; ] + [COZ] 
Kı[H2C0;] Kı K2[H2CO3] 


= [H2CO3] + 


m] mF 
Ki KK 
= icos(i + Ca ae) 
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Introducing the ionization fraction and a 


Ky KK 
ao=(1+ ta + Tas (6.22) 


we can express carbonic acid, bicarbonate and carbonate concentrations as a function 
of DIC, Kı, K2 and [H+]: 


{[H2CO3] = 2 (6.23a) 
Q 
DIC. K 
[c03] = orcs (6.23c) 
a . 


This ionization approach can not only be used for closed systems with known DIC, 
but also for open systems (without a priori information on DIC, but with known pH). 
In an open system, the carbonic acid concentration is set by Henry’s Law (Eq. 6.10), 
and via Eq. 6.23a one can then estimate DIC first and subsequently the bicarbonate 
and carbonate ion concentrations. 

Figure 6.1 shows the distribution of species contributing to DIC as a function of 
pH, a Bjerrum plot, calculated using the ionization fraction approach. Carbonic acid 
is the dominant species at pH values below the pK, bicarbonate dominates between 
the pK; and pK, values and the carbonate ion dominates at pH values above the 
pK2. Moreover, the pH dependent distributions of carbonic acid, bicarbonate and 
carbonate differ between freshwater and seawater because the many ion interactions, 
i.e., non-ideal behavior, in seawater shift equilibria to lower pH values. 

Case 2 Rainwater pH: an open system 

Carbon dioxide in rainwater is in equilibrium with the atmosphere, and we have 
five unknown aqueous species: [H+], [O0H7], [H2CO3], [HCO; | and [co;]. We 
can start the calculation using Eq. 6.10 from the known atmospheric partial pressure 
of carbon dioxide (420 jratm) and Henry’s Law constant (3.4 x 107? mol kg"! 
atm!) to arrive at a carbonic acid concentration of 1.42 x 1075 M (= 107485). 
Next, we rearrange the first equilibrium relation: 


_ [H*][Hco;] eee, 
Kı = CO] => [HCO; | ` [H ] = K, + [H2C0O;] 


= 107635 x 1074-85 = 10711-20 


and see that the product of [HCO; ] - [H+] = 107!!20, Next, we consider the charge 
balance (Eq. 6.15), 
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Fig. 6.1 Distribution of carbonic acid, bicarbonate and carbonate as a function of pH in ocean and 
freshwater for a total DIC of 2000 uM and a temperature of 15 °C 


[Ht] = [OH] + [HCO; ] + 2[CO; 


and realize that [HCO; | > [Coz] and [H+] > [OH™], because after the addi- 
tion of carbon dioxide pH should be lower than 7. The charge balance can thus be 
simplified to 


[H+] ~ [HCO;] ~ 10756 > pH = 5.6 


Pure rainwater (not impacted by atmospheric dust or pollution) is thus slightly 
acidic. Increasing atmospheric carbon dioxide levels to 1000 or 5000 matm would 
lower rainwater pH values to 5.4 and 5.06, respectively. 

Humans have released large quantities of SO, to the atmosphere due to burning 
of sulfur bearing fossil fuels for industrial activities and transport. In the air some of 
the SO, gas is oxidized to SO; and both these gases have dissolved in rain. This had 
led to the production of sulfurous (H2SO3) and sulfuric acid (H2SO.z). As a result, 
pH in rain declined to values down to 3, i.e., acid rain, in some industrial areas. Such 
low pH values have consequences for mineral equilibria, availability of essential 
nutrients, and organisms’ physiology. In the 1980/90s environmental measures such 
as sulfur removal from exhausts have been implemented resulting in a gradual return 
to more natural pH values, but some of the impacted ecosystems (e.g., forests) are 
still recovering from that perturbation. 
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Case 3 The pH of a solution in equilibrium with calcium carbonate and the 
atmosphere 

This case would correspond to a pool or shallow lake in a carbonate rock setting, 
or a coastal ocean with carbonate minerals at the seafloor. It also applies to travertine 
formation in rivers or stalagmite and stalactite formation in caves (see Box 4). The 
system is characterized by six unknown aqueous species: [H+], [OH™ ], [H2CO3], 
[HCO; | and [Coz] as in the previous case of rainwater, plus [Ca**] from the 
equilibrium with carbonate minerals. Consequently, we need six equations to solve 
the problem. The first four are the self-ionization of water (Eq. 6.6), the first and 
second acid-base equilibria (Eqs. 6.13 and 6.14) and Henry’s Law combined with the 
Pco, (Eq. 6.10). Since the water bodies are in contact and equilibrium with calcium 
carbonate minerals, the solubility product (Eq. 5.16) provides the fifth constraint: 


Ksp = [Ca] [C03] (6.24) 
The charge balance provides the sixth, final constraint. 
[H+] + 2[Ca**] = [OH] + [HCO; ] + 2[c03] (6.25) 


Notice that [Ca] and [Coz] are counted twice because of their double charge. 

With six unknown and six equations, the problem is well posed, but the mathe- 
matical solution is somewhat cumbersome. We will first present the full solution and 
then show that by making use of chemical insight, we can simplify the solution. 

Full solution: Since it is an open system in equilibrium with the atmosphere, we 
can start with Eq. 6.10: 


[H2CO3] = Ky - Pco, 


and insert it in Eqs. 6.19 and 6.21 to obtain expressions for bicarbonate and carbonate 
ions: 


Kı[H2C0;] _ KiKy- Pco, 


[HCO; ] = HJ A (6.26) 
2 K\K2[H2CO3] KıKıKų : Pco, 
[c 0; ]= [H+P = [H+P (6.27) 


Using Eqs. 6.24, 6.6, 6.26, and 6.27, we can rewrite the charge balance Eq. (6.25): 


Ks Ky Ki Kyi - Pco, Ki K2Ky- Pco 
T+ ness = TH] + m] +2 ny 2 (6.28) 


Multiplying by [H +F and re-arranging so that all terms are on the left-hand side 
results in 
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3 kalE T 
[=] t o - Peo K,,[H"] — Ki Ku: Pco, [H+] 
—2KıK:Kų - Pco, = 0 (6.29) 


This is a fourth-order polynomial in [H+] that can be solved by trial-and-error or 
numerical techniques. 

Approximate solution: For a system open to the atmosphere, the final pH is likely 
close to neutral (7 + 2), and we can thus simplify the charge balance (6.25) to: 


2[Ca**] = [HCO; ] (6.30) 
because protons and hydroxide ions are on the order of 1075 to 107° M and 


[HCO; | > [CO; ](see Fig. 6.1). Next, we rewrite the second acid equilibrium 
(Eq. 6.14) to isolate carbonate on the left-hand side: 


[H*][co;] 2-4 _ K2[HCO; ] 
= = 31 
= Eco] Hs T= e Ji 
Next, we substitute Eqs. 6.30 and 6.31 into the solubility Eq. 6.25 
K| HCO3 K: 
Ksp = [Ca (CO, |S Ksp = osco] Ecl = osea] Ta] 
(6.32) 


As above, equation (6.26) is used to express bicarbonate in terms of Pco, and 
[H+], and (6.32) is then rewritten as 


2 
ae Kə ven 


semos | 


Or after some further re-arrangements to isolate the proton concentration on the 
left-hand side: 


(6.34) 


[H+] = Poo? | SKA 


2K sp 


Solving this equation for a Pco, of 420 atm and using the standard values for 
PK, (6.35), pK2 (10.33), pK 4 (1.47) and the solubility product of calcite (pK = 8.48), 
yields a pH value of 8.18, consistent with observations for freshwaters in carbonate 
terrains. Moreover, it is also very similar to that of modern seawater despite seawater 
being a non-ideal solution. 

A more accurate way to calculate the pH for seawater in equilibrium with calcite 
and the modern atmosphere is to use equilibrium constants specific for seawater at 
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25 °C: pK” (5.87), pK3” (8.76), pK} (1.54) and the solubility product of calcite 
(p Ky,’ = 6.36). These seawater specific constants are shifted towards lower pK 
values, in other words higher K values (see Fig. 6.1). Modern seawater is estimated 
to have a pH of 8.09, consistent with observations. 

Atmospheric carbon dioxide levels have varied over the Earth’s history and 
Table 6.2 shows the pH value expected for freshwater and marine waters in equi- 
librium with calcite and the atmosphere at 25 °C for glacial, pre-industrial, modern, 
Eocene, and projected levels in 2100 for two global warming scenarios (2 °C and 4 
°C). 

Table 6.2 clearly shows that pH declines with increasing Pco, levels. This ocean 
acidification is also known as the other CO2 problem. The above calculations 
are based on full equilibration of seawater with both the atmosphere and calcium 
carbonate in the form of calcite. However, in the ocean, carbonate minerals are at 
the bottom of the ocean (depth of km’s) while gas exchange, and thus CO» inva- 
sion, occurs in surface waters. Ocean acidification for surface waters can better be 
calculated using the assumption of equilibrium gas exchange and constant alka- 
linity, i.e., no equilibration with calcite. In this case (last column), pH declines more 
steeply without this buffering effect of carbonate minerals, i.e., the surface ocean is 
more sensitive to atmospheric composition changes on timescales shorter than ocean 
turnover (needed to equilibrate with sedimentary carbonates). 

Case 4 The pH of a soda lake 

Evaporative lakes are usually alkaline and contain soda minerals such as natron 
(Na2CO0; - 10H20). The equilibrium reaction 


NayCO; - 10H20(s) <> 2Nat + CO} + 10H20 
constrains the carbonate ion and sodium concentration of the lake. Calculate the pH 
of the alkaline lake knowing that the dissolved sodium concentration ([Na*]) is 2.3 


g kg~!, that the lake water has a temperature of 25 °C and assuming equilibrium with 
the atmosphere, and a water density of 1 kg L~!. 


Table 6.2 Freshwater and ocean pH in equilibrium with calcite 


(Pco atm) | Freshwater pH with | Ocean pH with Ocean pH with TA 
calcite (Eq. 6.23) calcite (Eq. 6.23) | = 2750 uM 
Glacial 180 8.43 8.34 8.37 
Pre-industrial | 270 8.31 8.22 8.24 
Modern 420 8.18 8.09 8.09 
2°C-2100 475 8.14 8.05 8.05 
4°C-2100 800 8.00 7.90 7.86 
Eocene 1000 7.93 7.84 7.77 


Last column calculations with fixed TA rather than in equilibrium with calcite have been performed 
numerically 
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To solve this problem, we first identify the six unknowns species concentrations: 
[H+], [OH], [H2CO3], [HCO; ], [Coz | and [Nat]. Consequently, we need six 
equations to solve the problem. The first four are the self-ionization of water (Eq. 6.6), 
the first and second acid-base equilibria (Eqs. 6.13 and 6.14) and Henry’s Law 
combined with the Pco,(Eq. 6.10). The dissolved sodium concentration [Nat] is 
0.1 mol L~! and provides the fifth equation (mass balance). To close the system, we 
derive the charge balance: 


[H+] + [Nat] = [017] + [HCO; ] + 2[CO;"] (6.35) 


For an alkaline lake, [Nat | > [H+], and we can thus modify the charge balance 
and link it to the mass balance: 


[Na+] = [0H] + [HCO;] + 2[C03] = 0.1 (6.36) 


Since it is an open system in equilibrium with the atmosphere, we follow the 
approach for the previous case and use Eqs. 6.6, 6.26, and 6.27 to rewrite the charge 
balance (Eq. 6.36): 


Ky K, Ky - Pco, K,K2Kpx - Pco; 
+ ae) - 
[H+*] [H+] [H] 


=0.1 (6.37) 


and obtain an equation with all terms known except for [H+]. 
Substituting the appropriate values for the constants, 


10714 1071119 10721-53 
+ + 
[H+] [B+] [H+] 


and multiplying with [H +F we arrive at 
10-14. [H+] 41071129. [H+] 41072123 = 107 [H+] 
Which after re-arrangement turns out to be a quadratic equation 
[Ht] — 1971019. [H+] — 10-20-23 — 9 
and results ina [H +] of 1.1 x 107!° and thus a pH of 9.9. This estimate is approximate 


because of our neglect of non-ideal solution effects, i.e., using concentrations rather 
than activities. 
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Fig. 6.2 Conceptual figure of karst system with dissolution at surface and precipitation in caves 


Box 4 The Physical Chemistry of Karst 


Rainwater entering soils becomes enriched in carbon dioxide from respiration 
of the primary producers (trees and plants), and the microbial degradation of 
soil organic matter. These higher Pco, levels cause dissolution of carbonate 
bed rocks, and consequently dolines and sinkholes are formed, and streams 
and rivers may go underground. Soil waters will eventually recharge subsurface 
water reservoirs. When the groundwater drains into stream or rivers exposed to 
air, carbon dioxide will escape from the water to the air to re-establish equilib- 
rium. Consequently, the water becomes supersaturated with respect to calcium 
carbonates and travertine is formed. Similarly, when groundwater enters caves 
with lower Peo, levels, carbon dioxide will be transferred to the gas phase and 
stalagmites and stalactites may slowly precipitate (Fig. 6.2). 

Consider a typical soil Peo, level of 2000 ppm (because of extensive root and 
microbial respiration). Using equation (6.34) we can estimate that soil water 
in equilibrium with calcite (carbonate bedrock) would have a pH of 7.73. To 
calculate the dissolved calcium concentration in the soil, we must re-arrange 
Eq. (6.33) to isolate [Ca”*] on the left-hand side. 


(6.38) 


eo = Pana 


4K> 


Soil water recharging the groundwater would have a [ca] of 0.81 mM 
(32.5 mg L~!) and TA and DIC values of 1.63 mM each. If we assume that 
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[Ca], TA and DIC are transported conservatively, i.e., they do not change 
during transport, and that downstream cave/stream waters re-equilibrate with 
air with a Pco, level of 420 ppm, we can calculate the full carbon dioxide system 
in equilibrium. Cave and stream re-equilibrate with air and loose about 0.66 
mM DIC. The equilibrium pH, TA and [Ca°] would be 8.18, 0.98 mM and 0.48 
mM (19.3 mg L~!), respectively, because of calcium carbonate precipitation 
induced by the carbon dioxide transfer from water to air. 


Soil water (Peo, = 2000 ppm) | Cave/stream water (Pco, = 420 

ppm) 

pH 1.73 8.18 

[Ca] (mM) |0.81 0.48 

DIC (mM) 1.63 0.97 

TA (mM) 1.63 0.98 

[HCO;] (mM) | 1.62 0.97 

[co] UM) |41 6.8 
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Chapter 7 A) 
Redox Equilibria cere 


Abstract This chapter focuses on equilibria involving electron transfers (redox reac- 
tions). A systematic approach to identify electron transfers and to balance redox 
reactions is introduced. Redox potential and the Nernst equation are linked to the 
Gibbs free energy. 


Keywords Electron acceptor - Electron donor - Redox potential - Nernst 
equation - Oxidation state 


7.1 Introduction 


Redox reactions involve the transfer of electrons from an electron donor (reductant) 
to an electron acceptor (oxidant). The electron donor/reductant is then oxidized, 
while the electron acceptor (oxidant) is reduced. This ‘reductant is oxidized’ and 
‘oxidant is reduced’ terminology could be seen as confusing, and for clarity we 
will avoid the use of the terms reductant and oxidant and limit the use of the terms 
oxidation and reduction. Electron transfers are linked to energy transfers and provide 
the basis for most life on Earth. Heterotrophs, including humans, utilize organic 
matter (an electron donor) and combine it with an electron acceptor (oxygen) to 
respire and release the energy for maintenance and growth. Some microbes (e.g. 
chemolitho(auto)trophs) use minerals such as pyrite (FeS2) or reduced gases such 
as methane (CH4) as an electron donor to generate energy upon electron transfer to 
an acceptor (oxygen, metal oxides or sulfate). The exposure of reduced substances 
(organic matter, iron sulfide) to atmospheric oxygen (oxidative weathering) causes 
not only transfer of electrons (oxidation), but also release of protons. An extreme 
example is acid mine drainage: the oxidation of metal sulfides causes proton release 
with acidic conditions as a result. Pyrite (FeS2) is a major constituent of mine waste 
tailings and reacts with oxygen: 


4FeS) + 1502 + 10420 > 4FeOOH +8507 + 16H* 


© The Author(s) 2024 TT 
J. J. Middelburg, Thermodynamics and Equilibria in Earth System Sciences: An 

Introduction, SpringerBriefs in Earth System Sciences, 

https://doi.org/10.1007/978-3-03 1-53407-2_7 


78 7 Redox Equilibria 


This reaction generates 4 protons per mol pyrite oxidized and leads to acidification. 


7.2 Identifying Electron Transfer and Balancing Redox 
Reactions 


Redox reactions change the oxidation state of substances. An oxidation state, or 
oxidation number, is a hypothetical charge of an atom, if all its bonds to other atoms 
were fully ionic. It quantifies the excess or shortage of electrons. Oxidation numbers 
are integers because only discrete electrons can be transferred. The rules for assigning 
oxidation numbers are as follows: 


1. An atom in its elemental state has an oxidation number of 0: e.g., Na, O2, Ar 

2. An atom in a monoatomic ion has an oxidation number equal to the charge of 
the ion: e.g., Ca**, Fe?t, F7 

3. Oxygen has an oxidation state of —2, except when bonded to an oxygen or fluor 
atom. 

4. Hydrogen has an oxidation state of + 1, except when bonded as a hydride (—1) 
to metals. 

5. Alkaline and earth alkaline metals (columns | and 2 of periodic table) have an 
oxidation state of + | and + 2, respectively, and the halogens (F, Br, Cl, I) have 
an oxidation state of —1, except when bonded to oxygen or nitrogen. 

6. The sum of oxidation numbers is equal to the net charge for a polyatomic ion. 


For instance, in sulfuric acid, H.SO,4, sulfur has an oxidation number of + 6, 
while nitrogen oxidation numbers in nitrate (NO, ) and ammonium ( N Hy. ) are 
+ 5 and —3, respectively. Assigning oxidation states to sulfur is challenging as it 
may be —2, —1, 0, + 3, + 4 and + 6, depending on the molecular setting. It is 
important to recall that these oxidation numbers are hypothetical because one would 
infer that the charge of the two sulfur atoms in thiosulfate (S2 077) is + 2, while 
in fact one sulfur has a charge of —1 and the other + 5. Some common (Fe, Mn, 
Ti-bearing) minerals have elements with mixed oxidation number: e.g., in magnetite, 
Fe30;4, two of the iron atoms are + 3 and the other is + 2. 

Balancing reactions involving electron transfers requires a systematic approach 
involving two steps, first the identification of the electron donor and acceptor, hence 
the electrons transferred, and then the balancing of the elements. There are two 
methods: the half-reaction, or ion—electron method, in which the oxidation and reduc- 
tion processes are separated and each balanced before re-combining them again, and 
the oxidation state or whole-reaction method. The latter involves the following steps: 


(1) Assign oxidation state numbers to each of the atoms in the equation and write 
the numbers above the atom. 

(2) Identify the atoms that donate or accept electrons, i.e., that increase or decrease 
in oxidation number. 
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Identify the electron donating and electron accepting compounds on the left- 
hand side and quantify the number of electrons they donate/accept. 

(4) Balance the electron donating and accepting compounds on the left-hand side 
(5) Balance all elements, except O and H, on both sides of the equation 

(6) Balance oxygen on both sides of the equation through addition of H2O 

(7) Balance hydrogen on both sides of the equation through addition of Ht 


To illustrate this method, consider the following reaction that may occur in surface 
sediments or in groundwaters: 


FeS+ NO; > FeOOH+SO; + No 


First, we identify the oxidation numbers (oxygen and hydrogen are —2 and + 1, 
respectively) 


SE Da 5 +3 +6 0 

FeS+ NO, > FeOOH+ S0% + Np 
and recognize that Fe donates one electron to go from + 2 to + 3, S donates 8 
electrons to go from —2 to + 6 and N accepts 5 electrons to go from + 5 to 0. 
Accordingly, FeS donates 9 electrons (1 from Fe and 8 from S) and NO; accepts 
5 electrons (N from + 5 to 0). Balancing the electron donor and acceptor on the 
left-hand side thus implies that 5 FeS react with 9 NO, . 

5FeS+9NO, > FeOOH + SO; +M 

The rest of the balancing is simply a book-keeping, first, we balance Fe, S and N: 


5FeS + 9NO; > 5FeOOH +5850% +4.5Np 


Next, we balance oxygen atoms by adding three molecules of H20, in this case 
to the left-hand side: 


5FeS +9NO; +3H)O > 5FeOOH +5507 +4.5N> 


And finally, we balance hydrogen atoms by adding one proton, in this case to the 
right-hand side: 


5FeS + 9NO; +30 > 5FeOOH +5807 +4.5N, + Ht 
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7.3 Redox Potentials and the Nernst Equation 


Electron transfers are associated with energy transfer. Electrochemical cells (e.g., a 
battery) convert Gibbs free energy of a reaction into electrical energy and vice versa. 
Galvanic cells make use of spontaneous chemical reactions to generate electricity, 
while electrolysis requires electrical energy to drive a non-spontaneous reaction. 

Consider a strip of zinc metal in an aqueous solution of copper sulfate, the strip 
will become darker because of copper precipitates and the blue color solution will 
fade because of the spontaneous reaction consuming the colored Cu? 


Zn(s)+ Cu** (aq) > Zn* (aq) + Cu(s) (7.1) 


Zn donates two electrons which are accepted by Cu™ and the energy released 
goes into heat. If the electron donating (oxidation) and accepting (reduction) part 
of the reaction are physically separated but connected via a wire and salt bridge (an 
ionic solution allowing migration of cations and anions) to maintain charge balance, 
an electrical current will flow through the wire because of a difference in potential 
(E, expressed in V(olt)). Electrochemists split redox reactions in two half-reactions, 
or half-cells, one for the oxidation and one for the reduction process. For reaction 
7.1, these half-reactions are: 


Zn(s) > Zn?+ (aq) + 2e for the oxidation (electrons donating reaction) (7.2a) 
and 


Cut (aq) + 2e” — Cu(s)for the reduction (electrons accepting reaction). 
(7.2b) 


The potentials for these half reactions have been measured relative to that of 
protons/hydrogen gas conversions which has been set at 0 V for both reduction and 
oxidation: 


2H (aq) + 2e7 — Hb(g), reduction, electrons accepting reaction (7.3a) 


Hy(g) > 2H *(aq) +2e , oxidation, electrons donating reaction (7.3b) 


The standard reduction potentials E? (with electrons on the left-hand side and 
at SATP) have been tabulated and can be used to calculate the redox potential of 
reactions (Table 7.1). Oxidation potentials are the reverse of reduction potentials. 
For instance, the standard reduction potential for Zn and Cu are: 


Zn’ +2e° > Zn E? = —0.76 V 
(7.4a) 
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Cu?’ +. 2e7 > Cu E? = +0.34V (7.4b) 


The overall potential for the complete reaction 7.1 is then — (— 0.76) + 0.34 = 
1.1 V, resulting in a positive voltage as it should be for a spontaneous reaction. Redox 
potentials are intensive properties, i.e., independent of the amount, and additive. 

The standard redox potential E° is related to the standard Gibbs free energy via 
the Nernst expression: 


A,G° = -nF E? (7.5) 


where A,G° is the standard Gibbs free energy of a reaction, n is the number of 
electrons transferred and F is the Faraday constant (the electrical charge of 1 mol 
of electrons: 96,500 C (C) per mol e7). Recall from physics class that 1 J = 1 V-C. 
The minus sign in Eq. (7.5) is needed because spontaneous reactions have a negative 
Gibbs free energy and a positive redox potential. Accordingly, we can calculate the 
Gibbs free energy of reaction 7.1: 


A,G°? = —nF E? = —2 x 96500 * 1.1 = —212kJ 


Table 7.1 Standard reduction potentials at 25 °C 


Reduction Half-Reaction E? (V) 
Co** ag + eT > Co% (aq) 1.92 
H202 (aq) + 2 H* (aq) + 2 €7 > 2H20 1.77 
Cet (ag) + e7 > Ce (aq) 1.72 
Mn** (aq) +e7 —> Mn** (aq) 1.5 

Cr207?~ (aq) + 14 Ht (aq) + 6 e7 > 2 Cr (aq) + 7 H20 1.36 
O2(g) +4H* (aq) +47 > 2H20 1.23 
Ag (aq) +> > Ags) 0.80 
Fe™ (ag) + 7 —> Fe”* (aq) 0.77 
Cu** (ag) + 2 e7 > Cu í) 0.34 
Cu” (ag) + e7 — Cut (ag) 0.16 
2H* (aq) + 2&7 > H2 (g) 0.0 

Pb”* (ag) + 2 e7 > Pb) —0.13 
Sn” (ag) + 2 €7 > Sno —0.14 
Co* a +2e7 > Co (9) —0.28 
Cd?* aq) + 2 e7 > Cds) —0.40 
Fe% (aq) +2 €7 — Fe) —0.44 
Zn™ (aq) +2e > Zn (s) —0.76 
Mn™ (ag) + 2 €7 > Mn —1.17 
Mg” (aq) + 2e- > Mg -2.36 
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The standard redox potential is also related to the equilibrium constant because 
A,G° is linked to K: (Eq. 5.10): 


K — A, G° nFE? (7.6) 
= 6X = ex . 
P\ RT P\ RT 


Redox potentials, like Gibbs free energy changes, depend on temperature and the 
composition of the reaction mixture. For Gibbs free energy, we have (Eq. 5.9): 


A,G =A,G° + RTInQ 


where Q is the reaction quotient. Since A,G° = —nFE° and A,G = —nFE, we 
arrive at 


—nFE = —nFE° + RTInQ 


Dividing by —nF, this results in the Nernst equation 


Remo l hnos- logQ (7.7) 
~ nF 8 i 


2.303 RT 
nF 

Recall from math class that In 10 = 2.303 log 10 = 2.303. 

At 25 °C, the term ee has a value of 0.0592 V and thus 


.0592 
ee ja 
n 


gQ. (1.8) 


Returning to our initial reaction (7.1), we consider the case that [Zn] has a 
concentration of 0.5 M and [c u’t] = 0.01 M. The potential would then be 


E =E gQ =1.1 log = 1.05 V. 


2 [Cu] 


_ 9.0592 | 0.0592, [Zn] 
n 


The redox reaction is favorable. 

We conclude with articulating the internal consistency and interchangeability 
of equilibrium constants between different branches of chemistry. Equilibrium 
constants can be derived from: 

thermochemistry, i.e., AG$ data obtained on pure phases (Eq. 5.9): K = 


—AG? 
exp( RT ) 


Y pd 
the composition of solutions at equilibrium (Eq. 5.4) K = ca r 

eq Pe 
and redox potentials (Eq. 6): K = exp( uke” r 


In the absence of uncertainties, the equilibrium constants should be identical. 
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